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.  ABSTRACT 

>A  software  package  was  developed  to  perform  conceptual 
design  of  submarines,  using  the  Computervision  CGP-200X 
Designer  System*  a  turn-key  computer  aided  design 
hardware  and  graphics  software  system.  The  philosophy 
behind  the  software  package  is  to  keep  all  major  design 
decisions  under  the  control  of  the  design  engineer, 
rather  than  embedding  key  decisions  in  the  program 
algorithms. 

Modules  are  provided  for  calculating  weight  estimates, 
principal  character ist ics  and  envelope  geometry, 
resistance*  weight  and  moment  balance,  and  the 
equilibrium  polygon.  The  package  interfaces  with  a 
pressure  hull  design  module  developed  separately  in  an 
O.E.  thesis  by  Marvin  Meade.  Interactive  graphics  are 
used  where  appropriate. 

The  software  requires  a  knowledgeable  naval  architect  as 
the  user,  but  does  not  require  extensive  knowledge  of 
computers  or  computer  aided  design  systems.  ^ — —  -  — 
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INTRODUCTION 


GENERAL  DISCUSSION 


Computer  Aided  Design  (CAD)  has  become  an  important 
tool  m  many  industries  over  the  past  ten  years.  The 
capabilities  of  CAD  systems  are  growing  at  a  rate  which 
defies  efforts  to  catalog  them.  In  the  aerospace  and 
automotive  industries,  as  well  as  in  large  archi tectural 
firms,  CAD  has  become  the  standard  mode  of  design. 

In  the  area  of  naval  ship  design,  CAD  has  been 
adopted  more  slowly,  for  a  number  of  reasons.  The 
comparatively  low  production  rates  bring  the  issue  of 
cost  effectiveness  into  greater  question  for  ship  design 
agencies.  The  tremendous  complexity  and  scale  of  naval 
combatants;  makes  the  design  of  a  flexible  and 
comprehensive  package  of  software  a  formidable  task. 

Despite  these  problems,  a  number  of  packages  for 
conceptual  and  preliminary  design  of  surface  combatants 
have  been  implemented  in  both  government  and  private 


sector  design  organizations.  The  packages  used  by 
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gover'nm&nt  agencies  have  tended  to  be  batch-oriented 
ship  synthesis  programs  with  little  or  no  graphic 
output,  and  virtually  no  real-time  interactive  graphic 
design  capability  (the  notable  exception  to  this  is  the 
use  of  interactive  interior  ship  layout  software). 

This  lack  of  interactive  graphics  has  been 
partially  hardware  driven;  the  concern  for 
standardization  and  the  lengthy  process  of  government 
certification  and  procurement  cause  a  substantial  delay 
between  the  initial  availability  and  subsequent 
installation  of  new  hardware.  In  addition,  the  CAD 
process  is  foreign  to  many  of  the  prominent  designers  m 
the  very  tradi tion-or j tinted  world  of  naval  architecture. 

Many  design  procedures  in  naval  architecture  are 
somewhat  subjective  and  difficult  to  quantify. 

A  further  impediment  to  *he  use  of  CAD  has  been  the 
tendency  or  many  early  software  packages  to  usurp 
designer  perogative  by  having  key  decisions  in  the 
design  process  “hard-wired"  into  the  program  algorithms. 

Such  software  does  not  utilize  the  expertise  of  the 
designer,  and  frequently  precludes  trading  off  design 
parameters  to  achieve  optimal  designs. 

In  the  area  of  submarine  design,  the  above  problems 
are  further  accentuated.  The  volume  of  new  submarine 
designs  is  even  lower  than  that  of  surface  combatants. 

Accuracy  becomes  critical?  where  feet  and  tons  may  be 


acceptable  units  for  tolerance  in  large  surface  ship 
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designs,  inches  and  pounds  are  frequently  the  desired 
accuracy  level  in  submarine  design.  This  is  due  to  the 
hydrostatic  requirement  to  balance  both  surfaced  and 
submerged,  as  well  as  to  the  premium  attached  to 
internal  volume  and  deck  surface  area. 

Submarine  designers  are  fewer  m  number  than  their 
surface  counterparts,  and  they  often  have  many  design 
relat lonships  tied  to  subjective  criteria  (i.e. , 
“designer’s  eye")  developed  over  years  of  experience.  A 
program  which  does  not  allow  the  submarine  designer  to 
retain  this  subjective  design  flexibility  is  not  likely 
to  be  accepted  by  the  design  community. 


SUBMARINE  CADs  A  HISTORICAL  PERSPECTIVE 


Although  Computer  Aided  Design  has  not  been  a 
primary  tool  in  submarine  design,  numerous  programs  have 
been  developed  to  accomplish  analytical  chores  that  are 
subsets  of  the  overall  design  process.  These  include 
propel lor  design,  resistance  estimates,  hydrostatic 
calculations,  and  structural  design.  These  programs  have 
been  primarily  analytical  in  nature,  providing  no 
interactive  graphic  capability,  and  are  usually  written 
for  main-frame  computer  systems. 

Several  packages  have  been  written  to  achieve  a 
complete  iteration  of  the  design  process.  These  include 
the  CODESUB  program  l 1 3.  It  provides  numerical,  rather 
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than  graph  1 Ci  output.  CODESUB  was  developed  by  the 
Center  for  Naval  Analyses  in  order  to  aid  in  projecting 
future  submarine  design  character ist 1 cs  for  both  our 
navy  and  its  potential  adversaries.  This  program  was  not 
intended  to  be  a  detailed  tool  for  the  naval  architect? 
and?  unfortunately ?  has  many  key  design  parameters 
imbedded  in  the  source  code?  which  restricts  its  use  m 
trading  off  alternative  designs.  Other  pad  ages  "graft" 
together  sections  of  existing  designs?  matching  the 
largest  hull  diameter?  to  provide  a  composite  design. 


ASSET  C2D?  developed  by  Boeing  Computer  Services?  while 
primarily  a  surface  ship  design  package?  may  be  expanded 
to  include  a  submarine  design  module  in  the  future. 
ASSET  provides  some  low-level  graphics  output?  but  is 
not  a  full?  real-time  interactive  graphics  pad  age. 

Thus?  there  is  a  clear  opportunity  for  the 
development  of  a  submarine  design  package  incorporating 
the  use  of  the  full  capabilities  of  currently  available 
CAD  systems.  This  thesis  represents  an  exploration  of 
such  a  design  package. 


PHILOSOPHY  FOR  THE  THESIS  PROGRAM 


The  software  developed  in  this  thesis  accomplishes 
a  "first  pass"  through  the  submarine  design  process  at 
the  conceptual  design  level.  In  conjunction  with  the 
concurrent  Ocean  Engineer  thesis  written  by  Marvin 
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Meade,  the  package  includes  calculation  of  weight 
estimates,  envelope  and  pressure  hull  geometry, 
resistance  (speed  and  power),  basic  hydrostatics,  weight 
and  moment  balance,  and  determination  of  the  equilibrium 
polygon. 

The  primary  philosophy  of  this  package  is  to 
develop  it  for  a  user  who  is  a  competent  naval 
architect,  but  who  may  not  necessarily  be  an  expert  in 
the  use  of  computers.  Crucial  design  decisions  are  left 
to  the  user,  rather  than  being  embedded  in  the  program 
code.  Opportunities  are  provided  to  override  program 
algorithms  where  necessary  to  implement  the  desires  of 
the  individual  designer.  In  addition,  if  the  designer 
does  posess  some  knowledge  of  CAD  systems  or  computer 
programs,  it  is  easy  to  exit  the  package  at  appropriate 
points  and  tailor  the  design  to  provide  more  detail  or 
alternate  geometry  and  analytical  procedures. 
Interactive  graphics  are  used  to  provide  a  clear 
visualization  of  the  design  under  development.  All 
program  input  parameters  are  couched  in  standard  naval 
architecture  terms  to  facilitate  the  user’s  interaction 
with  the  package.  The  title  chosen  for  the  software 
package  is  CADSUB  (Computer  Aided  Design  of  Submarines). 

Input  to  each  program  module  is  explained  in  the 
chapter  documenting  the  particular  module.  Samples  of 
program  output  are  provided  in  Appendix  I  to  this 


thesis 


HARDUIARE  s 

THE  COMPUTERVISION  DESIGNER  SYSTEM 


GENERAL  DISCUSSION 


The  optimum  choice  of  hardware  for  implementing  the 
philosophy  of  this  thesis  was  one  of  the  several 
high-quality  "turnkey"  CAD  systems!  all  of  which  permit 
the  use  of  interactive  graphic  design  techniques. 
Graphics  packages  for  main-frame  systeiiis  are  available 
from  several  prominent  suppliers!  but  are  a  compromise 
at  besti  since  the  processors  they  run  on  are  not 
optimized  for  graphics-intensive  applications. 

Suitable  systems  are  marketed  by  several  suppliers! 


THE  COMPUTERVISION  DESIGNER  SYSTEM 


The  system  which  Computervision  provided  time  and 
instruction  on  is  their  basic  "Designer  System" ,  built 
around  the  CGP-200X  processor.  The  heart  of  the  system 
consists  of  the  processor,  a  high-speed  tape  drive,  and 
one  or  two  300  megabyte  hard  disk  units.  The  CPU  is 
available  with  several  increments  of  random  access 
memory  (RAM).  The  particular  system  used  for  this  thesis 
was  configured  for  approximately  1.3  megabytes  of  core 
memory.  The  architecture  of  the  CGP-200X  is  16-bit. 

The  Designer  System  utilizes  two  distinct  operating 
systems.  Basic  file  management  and  system-level 
housekeeping  chores  are  handled  at  the  "OS"  level,  which 
is  similar  in  use  to  the  operating  systems  on  typical 
mini-computer  installations.  For  the  graphics  operating 
system,  Computervision  uses  the  "CADDS  4-X"  environment, 
which  is  tailored  for  graphics,  and  automatically 
interfaces  with  a  powerful  data  base  management  system. 
Graphics  commands  are  entered  in  a  simple  verb-ncjn 
syntax,  such  as  "INSERT  LINE"  or  "INSERT  SPLINE*. 

Input/output  (I/O)  is  handled  by  several  devices. 
For  system  administration  and  text  file  input,  small 
alpha-numeric  terminals  are  provided.  For  graphics  I/O, 
high-resolution  color  or  monochrome  terminals  are 
provided  as  part  of  a  comprehensive  workstation,  which 


-12- 


includes  a  thermal  printer  for  working  quality  text  and 
graphics  printout.  Also  provided  at  each  workstation  are 
a  digitizing  tablet  and  a  display  control  device.  This 
device  controls  background  and  graphics  intensities?  the 
number  of  lines  of  text  displayed  <4  or  24  lines)?  and 
dynamic  control  of  the  display?  zooming?  rotation?  and 
scrolling.  Output  of  finished  quality  is  available  from 
line  printers?  color  pen  plotters?  and  black  and  white 
electrostatic  plotters.  Large  format  output  is  supported 
by  the  plotters. 

PROGRAMMING  LANGUAGES 


Several  languages  are  provided  for  use  with  the 
Designer  System.  The  primary  language  is  Fortran-S?  a 
subset  of  ANSI  Fortran.  System  subroutines  are  provided 
for  interacting  with  the  graphics  and  data  base 
management  operating  systems.  Programs  entered  in 
Fortran-S  execute  with  the  greatest  speed?  and  allow 
creation  of  new  graphics  commands.  Because  of  limited 
CPU  space?  the  utilities  required  for  this  level  of 
programming?  which  are  often  provided  on  main-frame 
systems?  must  be  separately  programmed  for  each  module. 
Such  utilities  include  linking  and  loading  capabilities. 
Because  of  the  length  of  training  required  to  do 
competent  programming  at  this  level?  this  thesis  was 


programmed  in  an  alternate  language. 


Two  macro  languages  are  provided  for  the  system. 
Both  allow  creation  of  execute  files  using  standard 
graphics  commands.  VARPR02  is  a  non-compi lable  macro 
language  which  supports  Fortran-like  I/O  and 
computational  statements.  It  is  a  flexible  and  effective 
languagei  but  executes  rather  slowly.  NEWVAR,  the  other 
macro  language  provided?  is  compi lable?  and  therefore 
executes  far  more  rapidly.  The  I/O  for  NEWVAR  is 
cumbersome  and  difficult  to  format.  All  I/O  must  be 


handled  in  text  string  form,  and  converted  within  the 
program  to  numerical  values.  Despite  this  shor tcoming? 
the  execution  speed  of  NEWVAR  led  to  its  adoption  for 
use  in  this  thesis.  Ultimately?  'The  algorithms  in  this 
thesis  should  be  translated  to  Fortran-S  for  maximum 
speed  and  I/O  flexibility. 


PRELIMINARY  UEIGHT  ESTIMATES 


IMPORTANCE  OF  HEIGHT  ESTIMATES 


The  algorithms  employed  in  this  thesis  assume  a 
weight-driven  submarine  design.  This  assumption  means 
that  the  envelope  size  calculated  for  the  submarine  will 
contain  enough  volume  to  displace  the  same  amount  of 
water*  when  submerged*  as  the  weight  of  the  envelope 
displacement  for  the  submarine.  If  the  desired  design  is 
known  or  suspected  to  be  volume-limited*  an  arbitrary 
excess  of  lead  may  be  specified  during  the  weight 
calculations*  thus  driving  the  program  to  calculate  a 
larger  envelope  size.  The  designer  may  translate  this 
additional  weight  into  volume  manually  to  allow  for  a 
volume-limited  design. 

Since  the  estimated  weight  data  will  drive  the 
basic  dimensions  of  the  submarine*  the  weight  estimating 
algorithms  must  be  made  as  accurate  as  possible.  The 
designer  must  carefully  consider  the  input  data  to  the 


weight  module*  and*  where  appropriate- 


may  decide  to 


override  the  programmed  algorithms  and  specify  the  value 


for  a  particular  weight  group.  The  program  is  designed 
to  automate  the  calculations?  and  not  the  design 
judgement?  for  the  user.  A  knowledgeable  user  may  easily 
modify  the  algorithms  themselves  to  suit  particular 


WEIGHT  ESTIMATING  ALGORITHMS 


The  accounting  system  chosen  for  weight  estimation 
is  the  SMBS  (Ship  Work  Breakdown  Structure)  convention 
used  by  the  Naval  Sea  Systems  Command  for  ship  weight 
records  C33.  This  system  consists  of  seven  major 
numbered  groups  for  all  ships.  For  submarine  weight 
reports?  several  special  weight  categories  are  added.  A 
brief  summary  of  the  weight  groups  is  as  follows: 

GROUP  DESCRIPTION 

I  Structural 

II  Propulsion 

III  Electrical 

IV  Combat  Systems 

V  Auxiliary  Systems 

VI  Outfit  and  Furnishings 

VII  Armament 

A1  Sum  of  I -VI I 

Lead  Total  lead 

A  A1  +  Lead 


Variable  Load 


(fuel?  stores?  etc.) 


■*  ^ ■  .*■>*'?** 


NS< 


Normal  Surface  Condition 


Reserve  Buoyancy  Main  Ballast 

Submerged  Disp.  Self  explanatory 

Free  Flood  Non-buoyant  flooded  volume 

Envelope  Disp.  Submerged  Disp.  +  Free  flood 

The  algorithms  chosen  for  the  weight  estimating 
module  take  standard  design  parameters  as  input*  where 
possible*  and  provide  an  output  in  terms  of  a  percentage 
of  A1  weights*  NSC*  or  submerged  displacement*  as 
appropriate.  These  percentages  are  then  combined  with 
the  discrete  weights  that  are  calculated  or  input  for 
selected  weight  groups  to  determine  the  Al  weight  total* 
lead,  NSC,  and  other  weight  parameters.  The  designer  has 
the  option  of  specifying  weights  for  groups  II,  IV,  VI, 
and  VII*  which  override  the  programmed  algorithms.  The 
following  section  discusses  the  chosen  algorithms  for 
calculating  each  weight  group.  Throughout  this  thesis* 
the  following  conventions  will  be  used  to  indicate 
mathmatical  operations* 

"+"  =  addition 

=  subtraction 
=  multiplication 
•/"  =  division 
"A"  =  exponentiation 

GROUP  Is 

Weight  group  I  is  determined  as  a  percentage  of 


NSC*  In  general* 

%Wl=/( depth,  material,  size). 

The  size  dependence  is  accomodated  by  calculating 
W1  as  a  percentage  of  NSC.  The  algorithm  implemented  in 
the  program  151 

7.W1=C1+C2*DEPTH. 

Cl  and  C2  are  dependent  on  hull  material  used,  and 
are  determined  parametr ical ly  from  data  on  past 
submarine  designs  C4,53.  The  DEPTH  parameter  is  a  user 
input  of  maximum  operating  depth  m  feet.  The  accepted 
design  practice  of  designing  for  a  collapse  depth  of 
150%  of  operating  depth  is  implemented  through  the 
chcice  of  Cl  and  C2. 


GROUP  II : 

Weight  group  II  is  calculated  directly,  in  tons, 
rather  than  as  a  percentage.  The  number  of  parameters 
involved  xn  calculating  W2  is  greater  than  for  any  other 
weight  group. 

W2=/C horsepower ,  propulsion  type,  battery  type  & 
capacity  ( capaci ty=endurance) 3 

The  horsepower  variable  for  a  nuclear  plant  is 
simply  shaft  horsepower,  while,  for  a  non-nuclear  plant, 
two  separate  horsepowers  are  needed?  shaft  horsepower 
(usually  the  horsepower  of  the  electric  final  drive 
motor)  and  charging  horsepower  (power  of  the  prime  mover 
used  to  drive  the  alternators  or  generators  to  provide 
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electrical  power  to  recharge  the  storage  batteries).  The 
weight  of  the  charging  prime  mover  is  a  strong  function 
of  technology  chosen?  ie  diesel*  wankel,  etc. 

The  calculation  of  storage  battery  weight  for  a 
non-nuclear  plant  is  dependent  on  battery  type*  capacity 
in  kilowatt-hours  at  a  high  discharge  rate  ("sprint")* 
and  capacity  in  KW-H  at  a  low  discharge  rate  (endurance 
power).  The  options  provided  for  battery  type  in  this 
program  are  nickel-cadmium,  improved  lead-acid  (German 
VARTA  type),  and  “standard"  lead  acid  (U.S.  Trident 
type).  Constants  determined  from  the  energy  densities  of 
these  battery  types  at  the  two  discharge  rates  are  used 
to  determine  battery  weight  (the  larger  of  sprint  or 
endurance  rate  battery  weights  is  chosen  as  the  final 
battery  weight). 

Electric  motor  weight  is  determined  from  a 
multiplier  based  on  a  modern,  ai* — cooled  DC  motor.  The 
option  exists  to  provide  other  multipliers  based  on 
liquid  cooling,  superconducting  technology,  etc. 

For  the  non-nuclear  plant,  an  additional  option 
exists  to  provide  a  lumped  adjustment  to  weight  group  II 
to  provide  for  unconventional  propulsion  technology 
(fuel  cells,  Stirling,  low-power  nuclear,  etc.).  This 
final  adjustment  allows  the  designer  maximum  flexibility 
to  tailor  group  II  without  making  modifications  to  the 
program  code. 

For  a  "straight"  nuclear  plant. 
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W2=<Cl*SHP>/[ log(SHP) 3AC2 
where  SHP=shaft  horsepower 

Cl  is  from  parametric  data. 

For  a  non-nuclear  plant, 

W2=WE+WM+W£+dW2 

where  WE=charging  prime  mover  weight 
UI1=propulsion  motor  weight 
UBsbattery  weight 
dW2=adjustment  to  W2  (user  input) 
WE=C3*CHP 

where  CHP=charging  horsepower 

C3  is  dependent  on  prime  mover  type. 
WM=C4*SHP. 

WB=MaxC ( C5*KWHS ) , ( C6*KWHE ) 3 
where  KWHS=capaci ty  for  "sprint" 

KWHE=capacity  for  endurance 

C5,C6  are  dependent  on  battery  type. 


GROUP  Ills 

Although  weight  group  III  would  intuitively  appear 
to  be  a  direct  function  of  installed  generating 
capacity,  a  study  of  historical  design  data  C5  3  reveals 
that  the  group  III  weights  are,  in  fact,  very  closely 
approximated  by  a  straight  percentage  of  A1  weights. 
This  anomaly  is  best  explained  by  the  large  amount  of 
this  weight  group  attributable  to  the  electrical 
distribution  system  spread  throughout  the  submarine. 


m 
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This  power  distribution  system  weight  is  a  direct 
function  of  the  size  of  the  submarine*  since  the  amounts 
of  cable*  connectors*  and  distribution  panels  vary  with 
size  of  the  vessel.  Weight  group  III,  consequently,  is 
calculated  as  a  percentage  of  A1  weights: 

*/.W3=Cl. 


GROUP  IV: 

Determination  of  group  IV  weights  is  very  difficult 
to  implement  in  an  “automatic"  mode.  The  variety  of 
combat  system  components  and  great  variation  in 
equipment  weights  makes  the  selection  of  an  algorithm 
extremely  difficult*  as  most  schemes  provide  accurate 
results  for  a  limited  number  of  combat  suites*  at  best. 


r  -•  V\ 
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GROUP  V: 

l.'eight  group  V  is  accurately  determined  as  a 
function  of  submarine  size,  since  the  auxiliary  systems 
required  for  a  vessel  are  directly  proportional  to 
displacement  and  internal  volume.  Group  V  is  thus 
calculated  as  a  percentage  of  A1  weights: 

•/.W5=Ci  i=l,2 

where  Ci  is  based  on  submarine  type  (attack  or 
ballistic  missile). 

GROUP  VI : 

Group  VI  weights  are  a  function  of  crew  size,  and, 
to  some  extent,  typical  mission  duration.  The  multiplier 
for  group  VI  weight  determination  is  based  on 
habitability  standards  commensurate  with  the  last  two 
classes  of  attack  and  ballistic  missile  submar ines, 
respectively.  The  weights  are  directly  calculated  in 
tons: 

W6=C1*NP 

where  NP=number  of  personnel  in  the  crew. 

GROUP  VII: 

The  estimation  of  armament  weights  is  nearly  as 
difficult  as  the  task  of  determining  group  IV  weights, 
since  the  possible  combinations  of  existing  and  future 
weapons  and  associated  launchers  are  infinite.  Again,  it 


is  highly  recommended  that  the  designer  specify  a  value 


for  weight  group  VII  based  on  his  off-line  determination 
of  the  weapons  payload.  As  was  provided  in  the  case  of 
combat  systems,  a  rough  estimate  of  group  VII  is 
available  from  a  percentage  of  A1  weights.  The 
percentage  constant  was  determined  from  historical  data, 
and  separate  values  are  provided  for  attack  and 
ballistic  missile  submarines. 

7.W7=C1. 

OTHER  WEIGHTS: 

The  amounts  of  lead,  variable  load,  reserve 
buoyancy,  and  free  flood  are  provided  as  fractions  by 
the  designer.  The  following  ranges  are  appropriate: 

Lead  (.07-.ll)*Al 

Variable  Load  < .04-.07)*NSC  (nuclear) 

(.08— . 19)*NSC  (non-nuclear) 
Reserve  Buoyancy  ( . 12-. 15 )*NSC 
Free  Flood  (. 05-. 1 )*( Submerged  Displ.) 

The  operator,  however,  may  specify  his  own  values 
as  desired. 

CALCULATION  OF  WEIGHT  SUMMARY 

When  either  percentages  or  actual  values  have  been 
calculated  or  specified  for  each  of  the  weight  groups, 
the  calculations  for  the  weight  summary  are  processed. 
At  a  minimum,  actual  weight  values  m  tons  have  been 
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determined  for  weight  groups  II  and  VI.  Additionally, 
the  designer  may  have  specified  weights  for  weight 
groups  IV  and  VII,  As  an  example  of  the  algorithms  used 
for  final  weight  calculations,  assume  that  only  W2  and 
W6  have  been  determined  in  tons.  The  other  permutations 
of  final  weight  calculations  are  similar,  but  allow  for 
dicrete  weights,  vice  percentages,  for  either  W4,  W7,  or 
both.  For  the  case  where  only  W2  and  W6  are  known,  the 
sequence  is  as  follows: 

(All  "7."  figures  used  are  converted  to  fractions? 
7.X=7.X/1B0> 

K 1  =  ( 7.W1+7.W  1  *7.LD )  /  ( 1  -7.VL )  +7.U3+7.W5 
where  7.LD=lead  fraction 

*/.VL=var  lable  load  fraction. 

A 1  =  ( W2+W6 )  /  <  1  -K 1  -7.W4-7.W7 ) 

W3=7.W3*A1 

W4=7W4*A1 

W5=7.w5*Al 

W7=7.W7*A1 

Lead=7.LD*Al 

A=A1+Lead 

NSC=A/<1-7.VL) 

Variable  Load=7.VL*NSC 
Main  Eal  last=y.RBY*NSC 

where  7RBY=fract ion  of  reserve  buoyancy 
Submerged  Disp.  =NSC+Mam  Ballast 
Free  Flood=’/CFF*Submerged  Disp. 
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where  7.FF=free  flood  fraction 

Envelope  Disp. =Submerged  Disp.+Free  Flood. 

The  calculated  weight  summary  is  displayed  on  the 
screen  with  a  summary  of  user  inputs?  at  which  time  the 
designer  may  accept  or  reject  the  results.  If  the 
results  are  rejected?  the  program  loops  back  to  the 
input  section  for  the  weight  module  and  another 
iteration  is  started.  If  the  results  are  accepted?  they 
are  written  into  an  output  file  for  record  purposes?  and 
various  input  and  calculated  parameters  are  written  into 
several  intermediate  files  for  use  in  later  program 
modules  as  the  design  is  completed. 

Because  of  the  high  degree  of  importance  attached 
to  the  calculated  weight  data?  future  refinement  of 
program  algorithms  should  include  a  priority  effort  to 
accomplish  two  major  tasks: 

(1)  Better  accuracy  of  calculated  weights. 

(2)  Further  subdivision  of  weights  (at  least 
partially  to  the  three-digit  SUBS  level). 
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PRINCIPAL  CHARACTERISTICS 
AND  ENVELOPE  GEOMETRY 


PRINCIPAL  CHARACTERISTICS 


The  results  from  the  weight  module  provide  the 
first  important  input  to  the  determination  of  the 
remainder  of  principal  character ist ics,  which  are 
geometric  m  nature.  Since  this  thesis  implements  a 
weight-driven  submarine  design  process,  the  next  design 
task  is  to  create  an  envelope  of  the  appropriate  volume 
to  contain  the  envelope  displacement: 

Envelope  Volumes=35*Envelope  Displacement. 

The  product  envelope  defined  by  this  module  of  the 
program  will  exhibit  the  following  geometric  principal 
character ist ics  (abbreviations  and  symbols  m 
parentheses  are  program  variable  names): 

Leng t h  Over a 1 1  ( LO A ) 

Length  of  Entrance  (Lf) 

Length  of  Run  (La) 

Mid  Body  Length  (Lm) 


Diameter  (D) 


»«■*** 
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Length/Diameter  Ratio  (L/D) 

Prismatic  Coefficient  (Cp) 

Bare  Hull  Surface  Area  (SF). 

In  addition,  once  envelope  geometry  is  defined,  the 
following  calculated  hydrostatic  character ist 1 cs  may  be 
found: 

Submerged  Longitudinal  Center  of  Buoyancy 
(LCBSUB) 

Draft  S  Normal  Surface  Displacement  (TNSC) 
Longitudinal  Center  of  Buiyancy  3  NSC  (LCENSC). 

These  parameters,  then,  when  coupled  with  the  input 
parameters  delineated  for  the  weight  estimation  module, 
constitute  a  set  of  principal  character ist ics  for  the 
design. 

DEFINING  THE  ENVELOPE 

A  number  of  algorithms  have  been  suggested  to 
calculate  the  geometry  of  the  envelope  for  submarines. 
These  range  from  a  "cut-and-paste"  approach  using  past 
designs  121  to  very  sophisticated  polynomial 
determination  of  offsets  C63.  In  between  these  extremes 
are  a  number  of  methods  utilizing  simple  thumbrules  and 
parametric  relationships  among  Cp,  D,  L/D,  and 
displacement  C4,7  3. 

Both  the  cut-and-paste  and  thumbrule  approaches  are 
somewhat  imprecise  m  allowing  the  designer  to  fully 


-27- 


control  the  geometry  of  the  resulting  hull  form.  The 
polynomial  methods,  while  very  precisely  controlled, 
require  the  designer  to  specify  design  parameters  which 
are  not  couched  in,  or  clearly  related  to,  traditional 
naval  architecture  vocabulary.  The  most  desirable 
algorithm,  then,  will  allow  strict  control  of  the 
resulting  geometry  while  accepting  input  parameters 
which  are  common  naval  archi tectural  terms. 

CHOSEN  ENVELOPE  GEOMETRY  ALGORITHM 

In  choosing  the  algorithm  implemented  in  this 
thesis,  two  restrictions  were  accepted  for  the  sake  of 
simplicity.  The  envelope  shape  is  confined  to  circular 
cross-sections,  and  an  analytically  smooth  envelope 
outline  is  calculated.  Figure  4.1  demonstrates  the 
nomenclature  and  subdivision  of  the  envelope  calculated 
with  the  chosen  algorithm. 

The  designer  is  prompted  for  the  following 

geometric  input  parameters: 

Length  of  Entrance 
Length  of  Run 

Forebody  Prismatic  Coefficient  CCpF) 

Aft  Body  Prismatic  Coefficient  (CpA) 

Maximum  Diameter  (D). 

The  algorithm  then  defines  an  envelope  consisting 
of  a  hemi-ellipt ical  forebody,  a  cylindrical  midbody. 


'Cs; 


*V  tf*.--1  '*"*•  •*3*  <***•  -_V  v-^-i  --, 


The  coefficients  Nf  and  Na  are  essentially  fullness 
coefficients  for  the  fore  and  aft  bodies*  respectively. 
In  order  to  facilitate  input  in  traditional  naval 
archi tectural  terms*  these  fullness  coefficients  are 
calculated  from  the  values  specified  by  the  designer  for 
prismatic  coefficients  of  the  fore  and  aft  body  sections 
of  the  envelope: 

Nf=  71 . 0477*CpFA5  +65. 8107*CpFa4  -483. 3037*CpFa3 
+587.4137*CpFA2  -281 . 7224+CpF  +49.5876 

Where  Nf  is  the  forward  fullness  coefficient. 

Na=*  379.  6546*CpAA5  -938. 4703*CpAA4 
+944.8853*CpAA3  -471 . 0872*CpAa2  +1 19. 1465*CpA  -11.3454 
Where  Na  is  the  aft  fullness  coefficient. 

These  polynomial  relationships  are  valid  for  values 
of  CpF  from  .5383-. 8720,  and  for  CpA  from  .3333-. 7111. 
These  ranges  substantially  exceed  those  for  actual 
submarine  designs.  The  relat ionships  were  derived  by 
using  known  (hand-verified)  equivalencies  of  prismatic 
and  fullness  coefficients  to  set  up  matrix  equations 
with  the  polynomial  coefficients  as  the  unknown  column 
vector.  Gaussian  elimination  was  then  used  to  solve  for 
the  unkown  coefficients.  Additional  hand-calculated 
equivalencies  were  then  utilized  to  verify  that*  within 
the  specified  ranges*  the  polynomial  relationships  yield 
results  correct  to  the  fourth  decimal  place. 

The  "X*  coordinate  for  the  overall  envelope  runs 
from  0  at  the  forward  perpendi cular  to  LOA  at  the  aft 


perpendicular-  In  the  formula  for  the  forebody 
hemi-el lipse,  a  local  X  coordinate  (Xf)  runs  from  0  at 
the  end  of  the  entrance  length  (Lf)  to  Lf  at  the  forward 
perpendicular : 


In  the  aft  body  parabolic  section,  a  local  X 
coordinate  is  again  used: 

Xa=  X— (Lf+Lm) 

Note  that  Lm,  the  mid  body  length*  is  a  calculated 
parameter,  and  not  an  input.  The  calculation  of  Lm  will 
be  described  later  in  this  chapter. 

The  hull  radius,  r(X),  for  the  elliptical  forebody 
is  calculated  as  follows: 

r(X)=  R*C  (lMXf/Lf  )ANf  )A(1/Nf  )  3 
where  R=  D/2 

The  hull  radius  within  the  mid  body  is  simply: 


The  hull  radius  within  the  aft  body  parabolic 
section  is: 

r(X)=  R*C 1  — ( Xa/La) ANa  3 

The  offset  r(X)  is  calculated  for  each  of  21 
stations  evenly  spaced  along  the  hull  length.  Obviously, 
r(X)  at  both  the  forward  and  aft  perpendiculars  is  zero. 
When  the  offsets  for  the  21  stations  have  been 
determined,  the  sectional  area  of  each  station  is 
calculated: 

SA  <  X ) =  Pi*r ( X  >  A2 


The  station  sectional  areas  are  then  integrated 
along  the  length  to  determine  the  volume  enclosed  by  the 


envelope: 

rLOA 

VOL=  /  SA ( X )  dX 

This  and  all  other  integrations  m  this  thesis  are 
done  numerically «  using  Simpson’s  Rule. 

On  the  first  pass  through  this  process,  the  mid 
body  length  is  zero,  and  LOA=  Lf+La.  The  resulting 
volume,  then,  is  the  total  volume  for  the  fore  and  aft 
body  sections  alone.  To  determine  the  proper  mid  body 
length,  the  volume  determined  above  is  subtracted  from 
the  required  envelope  volume,  and  Lm  is  calculated  as 
the  length  of  a  cylinder  of  radius  R  containing  the 
remaining  volume: 

DV=  Envelope  Volume- VOL 

Lm=  DV/(Pi*RA2> 

If  DV  is  negative,  an  error  message  is  printed,  and 
the  designer  is  prompted  for  new  geometric  input  values. 

After  Lm  is  calculated,  the  final  value  of  LOA  is 
determined: 

LOA=  Lf+Lm+La 

The  sectional  areas  for  the  resulting  hull  form  are 
calculated  at  the  new  21  stations,  after  which 
circumferences  and  bare  hull  surface  area  are 
determined: 

Cxrc(X)=  2*Pi*r  <  x) 


SF=  J  Circ(X)  dX 


The  overall  prismatic  coefficient  is  calculated  by 
comparing  the  envelope  volume  with  the  volume  of  a 
cylinder  (Vcyl)  of  length  LOA: 

Vcyl=  Pl«RA2*L0A 

Cp=  Envelope  Volume/Vcyl 

These  calculated  parameters,  when  taken  together 
with  the  input  values,  constitute  the  geometric 
definition  of  the  envelope. 


CALCULATION  OF  BASIC  HYDROSTATICS 


In  order  to  later  balance  the  submarine,  it  is 
important  to  accurately  determine  the  longitudinal 
centers  of  buoyancy  in  the  submerged  (LCBSUE)  and  normal 
surface  (LCBNSC)  conditions.  Calculation  of  LCBNSC 
requires  the  determination  of  the  draft  at  NSC.  LCB’s  in 
both  conditions  are  calculated  by  integrating  the  first 
longitudinal  moment  of  the  sectional  areas. 

Using  the  sectional  areas  for  the  fully-submerged 
hull  as  calculated  above,  LCBSUE  is  determined: 


( 


LOA 


LCBSUE-  /  X*SA(X)  dX  /VOLENV 


There  are  several  possible  approaches  to 
determining  the  draft  at  NSC.  When  the  calculations  are 
done  manually,  it  is  customary  to  choose  at  least  three 


points,  commonly  at  drafts  equal  to  D,  D/2,  and  D+R/2, 
and  determine  the  volumes  and  resulting  displacements 
for  each  of  these  drafts.  A  curve  of  displacement  versus 
draft  is  then  constructed,  after  which  the  draft 
corresponding  to  the  NSC  displacement  is  read  from  the 
curve.  While  this  method  is  expeditious  and  acceptable 
for  hand  calculat ions,  relatively  large  errors  may 
result  when  the  designer  is  not  experienced  enough  to 
draw  the  proper  curve  through  the  three  points.  Use  of 
additional  points  improves  the  accuracy,  but  the  manual 
calculations  are  tedious  and  time-consuming. 

With  the  speed  and  accuracy  of  the  computer 
available,  a  much  more  precise  determination  of  draft  at 
NSC  (TNSC)  is  possible.  A  potentially  important  factor 
that  is  often  overlooked  in  manual  methods  is  the  volume 
of  the  free  flood,  which  is  contained  in  the  envelope, 
and  may  be  a  signifigant  portion  (2-102)  of  the  envelope 
volume.  Study  of  existing  designs  indicates  that 
approximately  802  of  the  total  free  flood  volume  is 
submerged  at  TNSC.  Consequently,  the  "target" 
displacement  m  determining  TNSC  is  not  simply  NSC,  but 
rather  includes  this  portion  of  the  free  flood: 

NSCVOL=  35*(NSC+.8*Free  Flood) 

(for  802  free  flood).  Where 

NSCVOL  is  the  desired  volume  at  TNSC. 

The  method  chosen  to  determine  NSC  draft  is  a 
bracket  and  halve  iteration  on  the  variable  TNSC, 


beginning  at  TNSC=  3*R/2.  The  allowable  error  in  NSC 
volume  is  specified  as  0.52  in  the  current  program.  This 
tolerance  may  be  decreased,  with  additional  iterations 
required  for  convergence.  The  value  of  0.5X  s  chosen 
as  sufficient  for  the  conceptual  design  pht  .  while 
minimizing  the  number  of  iterations  required.  The 
algorithm  employed  is  as  follows: 

DV=  . 005*NSCVOL 
NPLUS=  NSCVOL+DV 
NMINUS=  NSCVOL-DV 
DR=  . 5*R 
TNSC=  r+dr 

For  each  trial  TNSC,  the  value  of  each  station 
radius,  r(X),  is  compared  to  DR: 

If  (DR)>r(X),  then  r(X>  and  SA(X)  remain 


unchanged. 


If  (DRXr(X),  then: 


SA ( X )  =  SA(X)  — C  r  ( X )  A2  *ARCOS(DR/r (X ) )  -DR 


*(r ( X ) A2  -DRA2) A.5  3 


(Zero  trim  is  assumed). 


ytOA 

VOL=  /  SA (X )  dX 


If  VOL>NMINUS  and  VOL<NPLUS,  then  the  correct  TNSC 
is  determined.  If  VOL<NMINUS,  then  DR  is  increased  by 
•5*DR.  Similarly,  if  VOL>NPLUS,  then  DR  is  decreased  by 
the  same  increment.  In  either  case,  the  algorithm  is 
repeated  until  the  value  of  VOL  is  acceptable.  When  the 


correct  TN5C  is  found,  LCBNSC  ic  calculated: 


LCBNSC=  /  X*SA(X)  dX  /NSCVOL 


The  calculated  results  are  displayed  on  the  screen 
for  the  designer’s  approval.  If  the  results  are  deemed 
unsatisfactory,  the  designer  may  reject  the  envelope 
geometry  and  return  to  the  input  section  for  another 
envelope  design  iteration. 

If  the  results  are  accepted,  all  of  the  specified 
input  values  and  calculated  data  are  written  to  an 
output  file  for  record  purposes,  and  appropriate  data 
for  other  program  modules  are  written  to  passing  files 
for  later  use. 


PLOTTING  THE  ENVELOPE 


An  accurate  plot  of  the  envelope  is  imperative, 
since  the  outline  displayed  on  the  monitor  will  guide 
the  designer  in  his  interactive  design  of  the  pressure 
hull.  The  21  stations  determined  thus  far  are 
insufficient  to  fit  an  accurate  curve,  since  the 
arbitrary  equal ly-spaced  stations  may  not  fall  at  the 
points  of  greatest  hull  curvature,  or  at  the  junctions 
of  the  mid  body  with  the  fore  and  aft  body  sections. 

To  adequately  portray  hull  curvature  and 
transitions  at  section  junctions,  the  program  calculates 
132  unequal ly-spaced  offsets  to  be  used  in  plotting  the 


envelope.  Offsets  are  most  closely  spaced  near  the 
forward  perpendicular  and  at  the  transition  points  of 
the  hull.  These  points  are  then  used  to  fit  a  Bezier 
spline  curve  which  forms  the  profile  outline  of  the 
envelope.  The  Plot  is  automatically  scaled  to  make  full 
use  of  available  screen  area,  and  both  vertical  and 
horizontal  scales  are  drawn  as  a  reference  for  the 
designer. 

The  two-dimensional  envelope  outline  is  rotated  to 
form  a  surface  of  revolution,  and  a  mesh  pattern  is 
superimposed  on  the  surface  to  aid  visualization  of  the 
resulting  three-dimensional  hull  form.  In  addition,  an 
isometric  view  is  displayed  below  the  profile  for 
perspective. 

The  resulting  display  of  two  views  of  the  outer 
hull  constitutes  the  graphic  output  of  the  envelope 
geometry  module.  The  designer  is  now  ready  to  proceed  to 
the  speed  and  power  calculations. 


SPEED  AND  POUIER 


GENERAL  DISCUSSION 


The  calculation  of  resistance  for  the  envelope 
geometry  developed  in  the  previous  module  is  a 
straightforward  analytic  process.  The  formulae  employed 
are  the  familiar  drag  estimation  relat lonsh ips  from 
basic  hydrodynamics.  A  proposal  for  future 
implementation  would  be  the  addition  of  a  propel lor 
design  module  to  accurately  predict  the  propulsive 
coefficient*  in  conjunction  with  the  envelope  geometry. 

To  facilitate  endurance  calculations  for 
non-nuclear  designs*  the  calculation  of  total  electrical 
load  at  each  speed  is  included.  This  total  load  figure 
must  include  all  non-propulsion  loads  (hotel  and  combat 
system)  as  well  as  those  required  to  propel  the 
submarine. 

With  the  data  passed  from  the  envelope  geometry 
module*  an  estimate  of  appendage  drag  may  be  made.  This* 
when  combined  with  hull  form  and  surface  area  data* 


completes  the  input  to  the  speed  and  power  algorithm. 


ESTIMATION  OF  APPENDAGE  DRAG 


Study  of  existing  submarines  reveals  that  appendage 
drag  comprises  thirty  percent  or  more  of  the  total  drag 
on  the  submarine  for  most  designs.  In  the  later  stages 
of  the  design  process,  detailed  consideration  should  be 
given  to  the  configuration  and  size  of  each  appendage  to 
improve  the  accuracy  of  resistance  estimates.  At  the 
stage  of  conceptual  design,  however,  details  concerning 
appendages  are  often  uncertain.  Many  designs  are 
formulated  without  final  selection  of  control  surfaces, 
for  example.  The  required  shape,  location  and  sizes  of 
surfaces  for  cruciform,  X,  and  inverted  Y  sterns  may 
exhibit  some  variance,  and  the  final  configuration  may 
not  be  chosen  until  much  later  in  the  design  process. 

In  view  of  this  uncertainty,  a  reliable  rough 
estimate  of  total  appendage  drag  is  desirable  for  the 
conceptual  design.  This  thesis  incorporates  the  method 
suggested  by  Bukalov  C8  3.  Bukalov  studied  numerous 
actual  submarine  designs,  and  developed  an  algorithm 
based  on  this  parametric  study.  This  algorithm  may  be 
reduced  to  a  single  formula  dependent  on  length  overall 
(LOA)  and  maximum  diameter  (D).  The  appendage  drag  <Da) 
is  then  calculated  as  follows: 


Da=  1 . 09065E-3*LO A*D+ 11.25 
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CHOSEN  ALGORITHM 


The  basic  relationships  implemented  in  this  program 
module  may  be  found  in  numerous  introductory  texts 
£9,103.  The  interactive  input  from  the  terminal  consists 
of  the  following  items: 

Propulsive  Coefficient  (PC) 

Non-propulsion  Loads  in  KW  (HL) 

The  remainder  of  the  required  input  is  passed  from 
the  envelope  geometry  module,  and  consists  of: 

Prismatic  Coefficient  (Cp) 

Length  Overall  (LOA) 

Maximum  Diameter  (D) 

Bare  Hull  Surface  Area  <SF) 

The  speed  and  power  data  is  calculated  and  stored 
in  an  array  indexed  on  speed  m  knots  (Vk)  from  zero  to 
forty  knots.  After  appendage  drag  is  estimated  in  the 
manner  previously  described,  the  calculations  proceed  in 
an  iterative  loop: 

Re  <  Vk  >  =  ( Vk* 1 . 689*L0A ) / 1 . 27908E-5 
Where  Re  is  the  Reynolds  Number. 

Cf=  .075/C log(Re)-23A2 

Where  Cf  is  the  frictional  drag  coefficient, 
calculated  by  the  1TTC  convention. 

Cr=  Cf *C 1 . 5* ( D/LOA ) A 1 . 5+7* ( D/LOA ) A3  3 
+. 002* ( Cp-. 6) 
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Where  Cr  is  the  residual*  or  form,  drag 
coefficient,  calculated  after  Jackson  [43,  and  the 
last  term  is  a  correction  for  submarines  with 
considerable  parallel  mid  body. 

Ct=  Cf+Cr+. 00025 

Where  Ct  is  the  total  drag  coefficient  for  the 
bare  hull,  and  the  final  term  is  the  correction 
factor  from  tow  tank  studies. 

HP ( Vk ) =  ( 1 /PC ) * . 00872*VkA3* ( Ct*SF+Da  > 

Where  HP(Vk)  is  the  drag  in  horsepower  for  the 
given  speed,  Vk. 

KW ( Vk ) =  . 7457*HP ( Vf  >  +HL 

Where  KW<Vk>  is  the  total  load  in  kilowatts, 
including  non-propulsion  loads,  for  the  speed,  Vk. 
These  calculated  results  for  drag  in  horsepower  and 
total  load  m  kilowatts  are  displayed  on  the  terminal 
screen  for  designer  inspection.  From  these  results,  the 
designer  may  determine  the  maximum  speed  of  the 
submarine  for  the  installed  shaft  horsepower,  as  well  as 
the  submerged  endurance  at  any  speed  for  installed 
energy  storage  capacity.  If  desired,  the  module  may  be 
re-run  for  'arious  values  of  propulsive  coefficient  and 
non-propulsion  loads,  until  a  satisfactory  result  is 
obtained. 

When  the  designer  accepts  the  calculated  values, 
the  results  are  written  into  an  output  file  for 
hard-copy  record  purposes. 
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PRELIMINARY  BALANCE 


GENERAL  DISCUSSION 


When  the  designer  has  completed  the  speed  and  power 
module,  he  will  proceed  to  the  design  of  the  pressure 
hull,  which  is  a  process  covered  separately  in  a 
concurrent  ocean  Engineer  thesis  by  Marvin  Meade  till. 
With  a  graphic  display  of  the  envelope  and  pressure  hull 
to  serve  as  a  visual  aid,  the  designer  will  enter  the 
subject  module  for  this  chapter  to  accomplish  a 
longitudinal  weight  and  buoyancy  balance. 

The  purpose  of  the  preliminary  balance  is  to 
determine  the  locations  of  lead  and  main  ballast  tanks, 
in  order  to  confirm  the  feasibility  of  the  pressure  hull 
design.  The  require-  i  *»ad  (Pb)  longitudinal  center  of 
gravity  (LOG)  is  «.  eved  by  balancing  the  longitudinal 
moments  of  all  weights  except  the  main  ballast  tanks  to 
create  an  attitude  of  zero  trim  in  the  normal  surface 
condition  (NSC). 

Poliowing  the  determination  of  lead  placement,  an 


additional  moment  balance  is  computed  to  locate  the  main 
ballast  tank  LCG  in  order  to  achieve  zero  trim  in  the 
submerged  condition. 

In  addition  to  the  longitudinal  balance,  a  vertical 
weight  balance  is  computed  to  provide  an  estimate  of 
transverse  stability  while  submerged.  This  estimate  is 
in  the  form  of  RG*  the  distance  from  the  vertical  center 
of  buoyancy  to  the  vertical  center  of  gravity  of  the 
submarine.  This  parameter  should  have  a  value  of  about 
one  foot  for  an  optimum  design.  To  simplify 
calculations,  no  trim  is  allowed  in  the  surfaced 
condition. 


CALCULATION  PROCESS 


The  first  set  of  input  data  required  for  balancing 
the  submarine  is  passed  from  previously  completed 
modules.  The  data  which  is  read  in  from  earlier 
calculations  is  as  follows: 

All  16  weight  categories  from  the  Weight  Module 


Longitudinal  Centers  of  Buoyancy  at  NSC  and 
submerged  (LCBNSC  and  LCRSUR) 

Maximum  Envelope  Radius  <R> 

Length  Overall  (LOA). 

The  first  set  of  keyboard  input  is  the  LCG  and  VCG 
information  for  the  seven  primary  weight  groups  and  the 
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variable  load: 


LCG( 1>-LC6(7) 


VCG< 1 )— VCG(7 ) 


Although  detailed  arrangements  have  not,  as  yet, 
been  addressed  in  the  design  process,  the  designer  must 
have  given  enough  thought  to  rough  arrangements  to 
provide  these  location  data  for  the  aggregate  weight 
subdivisions.  If,  when  detailed  arrangements  are 
completed  at  some  later  time,  weight  group  LOG  and  VCG 
locations  are  signi f igant ly  different  from  the  input 
provided  for  this  module,  the  balance  process  must  be 
repeated. 

The  prog.  *m  now  calculates  arrays  of  longitudinal 
and  vertical  moment  data: 

LM  < i ) =  LCG( i )*W( l ) 

VM<i)=  VCG(i)*W(i> 

Where  LM  and  VM  are  the  longitudinal  and  vertical 
moments  of  the  weight  groups,  respectively. 


LCGtLond.  A-l>  =  /  LMd)  /  /  W(i) 


VCG ( Con d.  A-l)=  /  VM(i)  /  >  W< i > 


The  designer  is  now  given  a  screen  display  of  the 
total  weight  of  lead,  and  is  asked  to  input  the  desired 
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amount  of  that  total  to  be  used  as  margin  lead  (Ml). 
This  lead  will  arbitrarily  be  located  at  LOA/2 
longitudinally  and  D/2  vertically.  The  remainder  of  the 
lead  is  considered  to  be  trim  lead  <TL).  The  required 
LCG  for  trim  lead  is  then  calculated: 

LCG(TL>=  CW(NSC)*LCBNSC  -LMlCond.  A-l )  -LM(VL) 
-ML*L0A/2  3/SL 

Similarly,  the  actual  VCG(Lead): 

VCG < Lead >=  <TL*6  +ML*R ) /W ( Lead ) 

Where  6  feet  above  the  keel  is  arbitrarily  chosen 
for  VCG(TL)  to  make  placement  in  the  circular  cross 
section  of  the  ballast  tanks  feasible- 
The  designer  is  then  asked  whether  the  required 
LCG C SL )  is  feasible.  If  the  response  is  affirmative,  the 


program  continues  on  to  main  ballast  tank  (MBT) 
calculations.  If  not,  the  designer  is  asked  for  new  LCG 
and  VCG  data  for  the  weight  groups  and  variable  load.  He 
is  prompted,  in  each  case,  with  the  previously  entered 
LCG  and  VCG  for  each  category.  Of  course,  if  reasonable 
adjustments  in  weight  group  locations  cannot  lead  to 
feasible  stability  lead  LCG  (feasible  being,  generally, 
between  the  physical  centers  of  the  mam  ballast  tanks), 
a  re-appraisal  of  the  pressure  hull  design  and  rough 
arrangements  is  necessary. 

When  the  lead  has  been  satisfactorily  located, 
calculation  of  required  MBT  LCG  is  accomplished 
(VCG(MBT)  is  arbitrarily  fixed  at  R): 


r 
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LCG( NSC )=  CLMCCond.  A-l >  +LM(Lead) 

+LM<  VL ) ]/W(NSC) 

LCG(MET>=  CW (Submerged )*LCBSUB 

-LCG ( NSC ) *W ( NSC ) D/W(MPT) 

The  designer  is  again  asked  to  decide  whether  the 
calculated  required  LCG  is  feasible,  with  similar 
results  in  program  flow  depending  on  the  answer. 
Feasibility  of  MET  location  is  dependent  on  the  desired 
geometry  of  physical  locations  for  MBTs  as  envisioned  by 
the  naval  architect.  The  continued  graphic  display  of 
envelope  and  pressure  hull  geometry  is  useful  in  aiding 
this  decision  process. 

The  program  now  sums  all  vertical  moments  for 
weight  groups,  lead,  variable  load  and  main  ballast  to 
arrive  at  an  overall  VCG  in  the  submerged  condition 
(VCGSUB).  An  estimate  of  submerged  stability  may  be 
calculated  and  displayed: 

BG=  R-VCGSUB 

Note  that  the  vertical  center  of  buoyancy  is 
assumed  to  be  at  R  (D/2)  for  the  submarine,  which  is 
basically  a  body  of  revolution.  Any  error  due  to 
buoyancy  of  appendages  will  be  small  and  conservative, 
since  off-center  buoyancy  items  (sail,  etc.)  are 
generally  above  the  geometric  center  of  the  envelope. 
Again,  the  designer  is  asked  to  accept  or  reject  the 
calculated  stability,  and,  again,  rejection  results  in 
looping  back  to  the  input  section  of  the  program  module. 


When  all  calculated  values  are  deemed  acceptable? 
the  input  data  and  results  are  written  to  an  output  file 
for  record  purposes,  and  the  locations  of  each  item’s 
center  of  gravity  are  plotted  on  the  pressure  hull 
graphic  display.  The  designer  is  now  ready  to  proceed  to 
the  calculation  of  the  equilibrium  polygon. 


THE  EQUILIBRIUM  POLYGON 


GENERAL  DISCUSSION 


One  of  the  most  important  sets  of  calculations  in 
the  submarine  design  process  is  the  determination  of  the 
equilibrium  polygon.  Any  submarine  experiences  changes 
in  weight  distribution  during  the  course  of  operations, 
as  items  in  the  variable  load  are  consumed.  A  system  of 
tanks  must  be  sized  and  located  in  such  a  uiay  that  this 
change  in  weight  distribution  and  amount  may  be 
compensated  to  allow  the  submarine  to  remain  in  proper 
balance  and  trim. 

A  detailed  list  of  the  variable  load  items  to  be 
considered  may  be  found  later  in  this  chapter.  In 
aggregate,  the  items  compose  major  groups  in  the 
categories  of  provisions,  ammunition,  stores,  fuel,  and 
other  consumable  fluids.  For  a  nuclear  or  other 
non-fossil-fuel  submarine,  the  boat  will  always  become 
lighter  as  items  are  consumed,  and  the  compensation 
consists  of  adding  appropriate  amounts  of  variable 


ballast  at  the  proper  locations.  In  the  case  of  a  diesel 
or  other  fossil-fuel  submarmei  a  paradox  exists  for  the 
designer.  As  fuel  is  consumed  from  the  fuel  ballast 
tanks  external  to  the  pressure  hull,  these  tanks  are 
compensated  with  seawater,  since  they  are  "soft"  tanks, 
and  rust  remain  full  of  fluid  to  maintain  a  zero 
differential  pressure  across  the  tank  structure.  The 
seawater  used  for  compensation  has  a  higher  density  than 
the  consumed  fuel  which  it  displaces.  Consequently,  the 
submarine  actually  becomes  heavier  as  fuel  is  consumed. 

The  tool  which  is  used  to  calculate  the  required 
locations  and  capacities  of  the  compensating  tanks  is 
called  the  equilibrium  polygon,  which  is  a  plot  of 
weight  versus  moment  (fore  and  aft)  of  the  weight. 
Figure  7.1  is  an  example  of  the  equilibrium  polygon  for 
a  submarine.  The  points  within  the  polygon  represent 
extreme  operating  conditions  for  the  submarine  in  terms 
of  the  weight  and  moment  required  to  balance  the 
conditions.  The  solid  lines  forming  the  "polygon"  are 
plots  of  the  weight  and  moment  added  by  progressive 
filling  and  emptying  of  the  major  compensating  tank 
groups.  These  groups,  which  may  be  comprised  of  several 
tanks  each  (the  aggregate  weight  and  moment  for  the 
group  is  plotted),  are  the  forward  trim,  auxiliary,  and 
after  trim  tank  groups.  If  all  of  the  operating  points 
for  the  submarine  are  contained  within  the  polygon,  the 
submarine  may  be  safely  compensated  for  all  expected 
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operational  conditions.  If  not*  the  designer  must  change 
the  location(s),  weight(s),  or  both*  for  the  appropriate 
group(s),  in  order  to  expand  the  polygon  to  include  all 
the  conditions. 

The  conditions  to  be  calculated  for  polygon  points 
are  specified  in  the  NAVSHIPS  Technical  Manual  (NSTM), 
Chapter  9290  C123.  The  extreme  conditions  are  more 
severe  than  an  actual  submarine  would  be  expected  to 
encounter*  and  represent  a  “worst  case-  situation  for 
weight  and  moment  balance.  These  extreme  conditions  may 
be  thought  of  as  representing  the  following  operational 
situations: 


HEAVY  #1:  Following  a  short*  fast  patrol,  during 
which  no  ordnance  it,  expended*  but  all  fuel  oil  is 


consumed. 


HEAVY  #2=  Same  as  the  previous  condition,  except 
that  only  the  fuel  oil  in  the  fuel  ballast  tanks  is 
consumed. 


LIGHT  #1  &  LIGHT  #2:  Two  variations  following  a 
short  patrol  with  no  fuel  consumed*  but  with  all 
ordnance  expended. 


HEAVY  FORWARD  #15  Following  a  patrol  where  only  aft 
ordnance  is  expended,  and  all  fuel  oil  is  consumed 


from  the  forward  fuel  ballast  tan! s 


HEAVY  FORWARD  #2:  Same  as  previous  condition? 
except  that  fuel  from  normal  fuel  oil  tanl s  is  also 
consumed. 

HEAVY  AFT:  Following  a  patrol  with  only  forward 
ordnance  expended?  and  with  fuel  oil  from  the  aft 
fuel  ballast  and  normal  fuel  oil  tanl'S  consumed. 

CONDITION  N:  The  "normal"  full-load  patrol 
condition  at  the  beginning  of  a  patrol.  For 
submarines  with  FBTs?  the  FBTs  a  specified  as 
fully  ballasted  with  seawater. 

CONDITION  M:  The  same  as  condition  N?  except  that 
submarines  with  FBTs  carry  a  full  fuel  load  m  the 
FBTs  (for  such  submarines?  this  condition  will  be 
lighter  than  condition  N). 

For  the  above  extreme  conditions?  water  density  is 
also  specified  within  typical  ocean  ranges?  so  as  to 
accentuate  the  worst  case  nature  of  the  operating 
conditions. 


s 
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THE  VARIABLE  LOAD 


The  items  that  comprise  the  category  referred  to  as 
the  variable  load  are  numerous  and  diverse  in  nature. 
Classification  of  the  items  is  again  aided  by  the 
guidelines  and  categories  specified  in  the  NAVSHIPS 
Technical  Manual. 

The  variable  load  is  subdivided  into  two  large 
sub-groups;  fixed  items,  not  expected  to  vary 
significantly  during  a  patrol,  and  the  truly  variable 
items,  routinely  expended  as  the  submarine  operates.  The 
•fixed"  portion  of  the  variable  load  consists  of  the 
following  items: 

CREW  AND  EFFECTS:  Actual  weight  of  the  submarine’s 
complement  of  personnel,  clothing,  and  associated 
personal  effects. 

SLEMs  OR  COMPENSATING  WATER:  The  weight,  for  a 
fleet  ballistic  missile  (FEM)  submarine,  of  the 
submarine  launched  ballistic  missiles  (SLEMs)  or 
the  compensating  water  in  their  absence.  SLBM 
compensating  water  exactly  matches  the  weight  and 
moment  of  the  missiles  themselves,  so  that  this 


category  may,  indeed,  be  considered  a  fixed  item. 


SANITARY  TANKS 2  The  weight  of  tanks  and  contents 


associated  with  waste  disposal. 

OXYGEN  CANDLES:  The  weight  of  these  devices,  used 
for  the  generation  of  oxygen  under  emergency 
conditions.  These  are  generally  a  very  small  weight 
item. 

LUEE  OIL  IN  SUMPS:  That  weight  of  lube  oil 
contained  in  the  oil  sumps  of  equipment. 

FIXED  CLEAN  FUEL  OIL:  That  weight  of  oil  carried  in 
non-compensated  tanks  which  are  maintained 
essentially  full  (i.e.,  shield  tanks  on  nuclear 
submarines) . 

As  can  be  seen  from  Reference  C12  3,  there  are  other 
items  specified  as  fixed  load,  but  they  are  minor  in 
nature,  and  beyond  the  level  of  detail  required  for 
conceptual  design.  For  example,  the  item  "depth  control 
tanks"  in  the  NSTM  is  not  included.  This  omission  will 
be  explained  in  the  discussion  of  fuel  load  subdivision 
m  the  following  section. 

The  items  considered  variable  in  nature  are  as 
follows: 

PROVISIONS  AND  STORES:  The  consumable  foodstuffs 


^ii 


(provisions)  ~.nd  spare  parts,  paper  supplies,  etc. 

(stores)  t.iat  the  submarine  crew  will  use  during  a  , 

patrol . 

REVITALIZATION  OXYGEN:  Oxygen  used  for  life  support 
(atmosphere  control)  during  lengthy  submerged 
periods.  On  newer  submarines  ( par t 1 cular ly 
nuclear),  a  method  of  o::ygen  generation  often 
obviates  the  need  for  oxygen  ban)  s. 

TORPEDOS,  MISSILES,  AND  AMMUNITION:  This  category 
is  se) f-explanatory.  All  expendable  ordnance  other 
than  SLBMs  is  included. 

WRT  TANKS:  This  item  thrs  "•nater  round  torpedo" 
tank  capacity,  used  to  compensate  for  the  presence 
or  absence  of  torpedos  in  the  tubes. 

RESERVE  ELECTROLYTE:  The  electrolyte  carried  to 
replenish  the  submarine  storage  battery. 

FRESH  WATER:  For  this  thesis,  this  is  a  composite 
category,  consisting  of  the  total  of  potable,  feed, 
and  battery  water. 

RESERVE  LUBE  OIL:  Lube  oil  carried  to  replenish  the 

i 

sump  lube  oil  during  the  patrol.  ? 

? 


FUEL:  Propulsion  fuel  for  the  submarine-  See  the 
discussion  that  follows  for  the  sub-categories  of 
fuel  and  their  significance  to  the  polygon 
computations. 

Fossil^fuel  submarines  present  a  variable  lead 
parodox  for  the  designer-  Fuel  which  is  carried  external 
to  the  pressure  hull  (the  majority  of  the  fuel  load  for 
most  "conventional**  subr  arines)  must  be  replaced  by 
seawater  as  it  is  burned,  since  these  external  "fuel 
ballast  tanks’*  (FBTs)  are  not  constructed  to  withstand  a 
high  pressure  differential  (i.e.,  they  are  "soft 
tanks").  The  seawater  which  is  used  to  compensate  these 
tanks  is  heavier  than  the  fuel  it  replaces-  Thus,  as  the 
submarine  consumes  fuel,  it  gets  heavier,  rather  than 
1 ighter • 

To  provide  a  margin  of  safety  at  the  base  of  the 
equilibrium  polygon,  some  means  must  be  provided  to 
remove  weight  as  fuel  is  burned.  In  earlier  designs,  a 
"safety  tank"  was  provided,  of  sufficient  capacity  that 
when  it  was  emptied,  enough  water  weight  was  pumped 
overboard  to  compensate  for  the  added  seawater  in  the 
FBTs  as  fuel  was  consumed.  The  safety  tank,  then,  was 
full  at  the  beginning  of  the  patrol,  and  would  be  empty, 
or  nearly  so,  when  all  fuel  had  been  removed  from  the 
FBTs-  For  purposes  of  this  thesis,  the  safety  tank 


concept  was  not  used,  since  the  weight  and  volume  of  the 
tank  and  it’s  contents  is  a  non-productive  load  item. 
However,  a  means  to  compensate  for  the  FET  weight 
addition  is  necessary. 

The  method  chosen  for  implementation  in  this  thesis 
is  the  provision  for  a  variable  fuel  oil  (VFO)  tank.  The 
VFO  tank  is  internal  to  the  pressure  hull,  and  does  not 
require  seawater  compensation  as  fuel  is  consumed.  In 
actual  practice,  this  tank  may  be  ballasted  and  treated 
as  a  dual-purpose  tank,  acting  as  an  additional  trim  or 
auxiliary  tank.  For  computational  simplicity  within  the 
program  module,  however,  the  VFO  tank  is  treated  as  an 
uncompensated  clean  fuel  oil  tank.  An  operational  mode 
of  employment  is  assumed  whereby  the  fuel  in  the  VFO 
tank  is  consumed  first  during  the  patrol,  before  the 
FBTs  are  utilized.  Thus,  the  size  of  the  VFO  tank  must, 
at  a  minimum,  be  sufficient  to  compensate  for  the 
difference  m  fuel  and  seawater  densities  for  the  FET 
capacity.  This  dictates  a  minimum  capacity  of  23%  of  the 
total  fuel  load.  The  designer  is  asked  to  input  the 
percentage  of  total  fuel  to  be  allocated  to  the  VFO  tank 
(or  tank  group,  as  arrangements  may  dictate).  If 
sufficient  internal  volume  is  available,  it  may  be 
advantageous  to  use  a  figure  greater  than  23%,  in  order 
to  provide  an  adequate  margin  at  the  base  of  the 
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CALCULATING  THE  POLYGON 


The  program  module  execution  begins  by  reading  the 
weight  group  values*  submerged  LCB*  and  LCG  of  the  main 
ballast  tanks  from  data  files  generated  by  the  weight 
estimating  and  balance  modules.  The  total  load  to 
submerge  the  submarine  is  calculated* 

Load  to  Submerge^  Submerged  Displacement-  Condition 
A  Weights. 

The  total  variable  load  weight  is  displayed*  and 
the  designer  is  prompted,  item  by  item*  for  the  weight 
and  LCG  of  each  variable  load  item.  As  each  weight  and 
LCG  are  entered,  the  balance  of  the  wei9ht  in  the 
variable  load  account  is  displayed.  Items  which  are 
commonly  distributed  in  two  sub-items*  fore  and  aft*  are 
input  as  two  separate  weights  and  LCGs.  Examples  of  such 
items  are  torpedos*  missiles*  and  fuel.  The  designer 
should  have  given  some  thought  prior  to  this  input 
sequence  to  the  subdivision  of  the  variable  load.  The 
envelope  and  pressure  hull  are  displayed  at  the  top  of 
the  screen  during  this  stage  of  the  mcdule*  to 
facilitate  estimating  the  LCG  for  each  item. 

The  last  item  to  be  entered  is  fuel.  The  balance  of 
the  variable  load  account  at  this  point  is  assumed  to  be 
all  fuel  (obviously*  the  balance  should  be  zero  for  a 
non-fossil-fuel  submarine).  The  total  fuel  load  is 
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displayed,  and  the  designer  is  asked  to  enter  the 
fraction  to  be  allocated  to  the  VFO  tank.  He  is  then 
asked  to  provide  LCGs  for  the  forward  and  aft  FETs  and 
the  VFO  tank. 

With  the  input  of  all  items  completed,  the  program 
displays  a  variable  load  summary  of  all  item  weights  and 
LCGs.  If  the  summary  is  satisfactory,  the  program 
proceeds  with  calculations;  if  not,  the  designer  is 
returned  to  the  input  sequence,  where  the  prompt  for 
each  item  now  includes  a  display  of  its  current  value. 
If  this  value  is  to  remain  unchanged,  the  operator 
simply  hits  the  "RETURN"  key;  otherwise,  a  new  value  may 
be  specified.  This  iterative  input  procedure  is  repeated 
until  the  designer  is  satisfied  with  the  variable  load 
summary. 

When  the  input  of  all  variable  load  items  is 
complete  and  accepted  by  the  operator,  the  program 
calculates  the  water  to  balance,  and  the  associated 
moment  of  that  water,  for  each  of  the  nine  specified 
polygon  points.  Note  that  conditions  M  and  N  will  be 
identical  for  nuclear  submarines.  Condition  M  is  only 
defined  :or  submarines  with  FBTs. 

A  computation  of  one  of  the  equilibrium  conditions 
would  proceed  in  the  following  manner.  Prior  to 
computing  the  individual  conditions,  the  arms  and 
moments  for  each  of  the  variable  load  items  are 


determined 


Item  Arm=  Submerged  LCE  —Item  LCG 

Item  Moment =  Item  Weight  * I  tern  Arm 

In  addition*  factors  of  proportional ity  are 
determined  for  light  and  heavy  seawater.  For  each 
condition,  the  sum  of  all  variable  load  items,  n  the 
fractional  amount  specified  by  the  NAVSHIPS  Technical 
Manual,  is  computed.  Similarly,  the  total  moment  of  the 
specified  items  is  determined.  The  required  water  to 
balance  and  moment  to  balance  may  now  be  calculated: 

Water  to  Ealance=  Load  to  Submerge  -MET  -Sum  of 

Item  Weights 

where  MET=  Main  Ballast  Weight 

Moment  to  Balance-  -MET  Moment  -Sum  of  Item  Moments 

If  required,  the  water  and  moment  to  balance  are 
multiplied  by  the  proportionality  factor  for  light  or 
heavy  seawater.  Each  of  the  nine  conditions  is 
calculated  in  a  similar  fashion,  resulting  in  arrays  of 
weight  values  in  tons  and  moments  in  foot-tons  (with 
negative  moment  values  indicating  an  aft  moment). 

Several  "housekeeping"  chores  are  performed  at  this 
point  in  the  module.  The  display  of  envelope  and 
pressure  hull,  and  associated  data  base  attributes,  are 
filed  under  the  part  name  MIT.HULLOUT,  for  later  recall 
to  plot  desired  views  of  this  graphic  screen.  The 
maximum  values  of  water  and  moment  to  balance  are 
examined  to  determine  the  required  scale  for  the  polygon 
plot.  A  new  part  is  activated  under  the  name  MIT.POLYOUT 


for  graphically  displaying  the  polygon.  Depending  on  the 
scale  required)  one  of  several  pre-constructed  "form" 
drawings  is  activated,  providing  the  axes  and 
appropriate  scale  information  for  the  polygon  under 
consideration. 

The  program  now  inserts  and  labels  each  of  the 
equilibrium  condition  points  on  the  graphic  display. 
This  gives  the  designer  a  clear  display  of  the  limits 
required  for  his  polygon  design.  An  important 
considerat ion  at  this  point  is  whether  an  adequate 
margin  remains  at  the  base  of  the  polygon.  Generally,  if 
the  lowest  weight  value  is  less  than  15  tons  from  the 
abscissa,  the  designer  should  seriously  consider 
re-structuring  the  variable  load,  as  necessary,  to 
achieve  this  margin.  This  may  be  accomplished  by 
returning  to  the  input  section  of  the  module.  An  example 
of  such  a  change  would  be  to  place  more  of  the  fuel 
load,  volume  permitting,  in  the  VFO  tank,  for  a 
fossil-fuel  submarine. 

The  operator  is  now  asked  to  provide  the  LCGs  of 
the  Forward  Trim,  Auxiliary,  and  After  Trim  tank  groups. 
He  is  then  prompted  for  the  capacities,  in  tons,  for 
each  of  the  tank  groups.  The  program  determines  the  arms 
and  moments  of  each  group  m  the  same  fashion  previously 
delineated  for  the  variable  load  items. 

Lines  are  now  inserted  on  the  display  to  represent 
the  successive  filling  of  Forward  Trim,  Auxiliary,  and 


After  Trim  tanks,  fol lowed  by  the  successive  draining  of 
the  tanks  in  the  same  order.  These  lines  form  the 
polygon,  and  all  equilibrium  points  must  fit  within  its 
boundaries  if  the  tank  locations  and  capacities  are 
adequate.  Obviously,  a  polygon  whose  boundaries 
substantially  exceed  the  extent  of  the  equilibrium 
points,  while  safe,  represents  a  waste  of  internal 
volume,  and  should  not  be  accepted. 

The  designer  now  has  the  option  to  accept  or  reject 
the  polygon,  as  plotted.  If  it  is  rejected,  the  polygon 
boundaries  are  erased,  and  the  designer  may  specify  new 
tank  group  capacities,  LCGs,  or  both,  after  which  a  new 
polygon  is  plotted.  This  Process  may  be  repeated  as  many 
times  as  are  required  to  achieve  an  acceptable  polygon. 

When  the  polygon  results  are  satisfactory,  the 
polygon  is  plotted  on  the  electrostatic  plotter,  the 
part  MIT. POLYOUT  is  deleted,  and  the  variable  load 
summary,  tank  locations,  and  tank  capacities  are  written 
to  an  output  file  for  a  permanent  record  of  the  final 
input  data  set.  One  iteration  of  the  submarine 
conceptual  design  process,  as  supported  by  this  thesis, 
is  complete. 


SUMMARY  AND 


RECOMMENDATIONS  FOR  FURTHER  RESEARCH 


This  thesis  constitutes  a  first  exploration  of  the 
use  of  interactive  graphics  for  conceptual  design  of 
submarines.  The  resulting  product  is  useful,  flexible, 
and  accurate  for  conceptual  design  activity.  It  is  by  no 
means  a  finished  product  in  the  commercial  sense,  but 
provides  a  tool  to  build  on  and  to  use  confidently  for 
such  applications  as  student  design  projects. 

The  power  of  real-time  graphic  display  of  the 
design  during  its  creation  is  a  major  factor  in  the 
worth  and  potential  of  the  program.  Many  hours  of 
laborious  manual  drafting  are  avoided  by  the 
high-resolution  graphic  output.  In  addition,  a  user  with 
even  a  rudimentary  knowledge  of  the  Computervision 
system  can  go  far  beyond  the  capabilities  of  the  basic 
software  package,  thus  greatly  enhancing  its  value  as  a 
design  tool. 
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The  program  is  highly  interactive;  the  user 
controls  all  important  design  decisions.  The  designer 
must  have  a  solid  knowledge  of  naval  archi tecture,  which 
is  appropriate  for  any  ship  design  program.  Subjective 
design  judgements  must  remain  in  the  province  of  the 
user  for  a  credible  result,  and  that  philosophy  has  been 
rigorously  adhered  to. 


RECOMMENDATIONS  FOR  FURTHER  RESEARCH 


There  are  nearly  unlimited  opportunities  for 
continued  work  on  this  project,  and  in  this  general 
area.  This  thesis  has  scratched  the  surface  m  an  area 
ripe  for  research  and  development. 

The  first  obvious  areas  encompass  those  portions  of 
the  design  process  which  were  excluded  from  the  current 
project.  These  include  structures,  propeller  design  and 
detailed  interior  arrangements. 

Submarine  structural  design  could  benefit  greatly 
from  treatment  in  the  CAD  environment.  Visualization  of 
complex  structural  relationships  is  intrinsically 
clearer  with  the  use  of  real-time  graphic  display.  In 
addition,  the  use  of  powerful  finite  element  packages, 
with  automatic  mesh  generation,  is  possible  on  many  CAD 
systems,  including  Computervision.  These  packages 
display  stress  levels  m  color,  and  will  magnify 
displacements  on  command  for  rapid  visual  analysis. 
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Arrangement  of  interior  details  is  feasible  on 
these  systems,  including  automatic  interference  checking 
and  layout  of  piping  and  electrical  distribution 
networks. 

Propeller  design  may  be  facilitated  by  graphic 
display  of  the  flow  regime  during  analysis.  Much 
analytical  work  remains  in  this  area,  and  research  into 
links  with  mam-frame  systems  and  large  data  bases  is 
also  needed. 

In  addition,  the  translation  of  this  package  into 
the  Fortran-S  language,  with  attendant  programming  of 
linking  and  loading  files,  would  greatly  improve 
execution  speed  and  I/O  flexibility.  The  creation  of 
customized  graphics  commands  would  add  more  power  and 
enhance  the  ergonomics  of  the  package. 

Further  work  m  the  subject  areas  of  weight 
estimation,  hydrodynamics,  and  animated  display  of 
control  system  performance  is  required  to  derive  the 
maximum  benefit  from  submarine  CAD.  All  of  these 
subjects  have  been  separately  addressed  in  other 
software  packages,  but  would  be  valuable  additions  to  an 
integrated  design  program. 

With  increased  interest  and  design  volume  in  the 
submarine  design  field,  research  funding  is  both 
necessary  and  clearly  justified  by  the  potential 
improvements  to  the  speed  and  accuracy  of  the  design 
process.  The  overall  area  of  computer  aided  design  of 
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submarines  can  support  major  research  projects  m  the 
futurei  benefitting  both  academic  and  commercial  design 


activities. 
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MIT. AB CD. HEIGHT 
3-28-64  13*37*44  FUTIL  6.18 


UNEIGHT  ESTIMATE  FOR  13.461  DESIGN  PROJECT* 
2! 

3* TYPE*  ATTACK 

4!PRoPULSlON*  MODIFIED  HANKEL  ELECThIC 
51BATTERY  TYPE*  NICKEL  CADMIUM 
61SPRINT  KH-H*  3300 
7! ENDURANCE  KH»H*  8000 
8JBATTERY  HEIGHT*  165.631 
0! ENGINE  KnRSEPoHErt*  4500 
10 'ENGINE  A  ASS'"*;.  EQUIP.  HEIGHT*  22.86 
USMOTOR  &  AS^OC.  EQUIP.  HEIGHT*  32.25 
l2!DISCRFTc  ADDITION  To  GROUP  2  HEIGHT*  80 
13!MAX.  OPERATING  DEPTH  IN  FEET*  700 
14* 

I* :LEAD  FRACTION*  0.1 

*6! VARIABLE  LOAD  FRACTION*  0.080 

17! RESERVE  BUOYANCY  FRACTION*  0.J25 

18! FREE  FLOOD  FRACTION*  0.08 

tO! 

20! 

21 '.CATEGORY  HEIGHT 
22! 

23! GROUP  1  403.318 

24*CR0UP  2  300.741 

25JCR0UP  3  17.4418 

26! GROUP  4  05 

271GROUP  5  137.043 

28! GROUP  6  117.3 

20! GROUP  7  85 

30JC0ND.  A I  1245.84 

311LEAD  124.584 

32JC0ND.  A  1370.43 

33! VAR.  LOAD  133.884 

34!N$C  4504.31 

35! MET  188.030 

36!SUB.  DISP.  1692.35 
3 7 '.FREE  FLO.  135.388 
381ENV.  DISP.  1827.74 


f 
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M1T.A8CD.GE0M0UT 

3-28-B4  13*38*56  FUT1L  6.18 


1  JGEOMETRY  OUTPUT  FOR  13.461  DESIGN  PROJECT* 
2! 

3! 

4JLENGTH  OVERALL*  130.947 
5JMID-B00Y  LENGTH*  74.0465 
6!OVERALL  PRISMATIC*  0.708368 
7! LENGTH/DIAMETER*  5.1832 
8!SUBMERGED  LC8  I  AFT  OF  FP)*  66.2480 
O'DRAFT  •  NORMAL  SURF.  COND.»  21.0585 
10'LCB  •  NORMAL  SURF.  COND.  *  66.3320 
II! 

12! RECAP  OF  INPUT  DATA* 

J3! 

14! ENVELOPE  DISPLACEMENT  *1827.74 

I5!£NTRANCE*  25 

16! RUN  *  40 

I7!DIAMET£R*  27 

13‘FND  PRISMATIC*  0.6 

IO!AFT  PRISMATIC*  0.55 


MlT.&BCD.OFFOUT 
3-28-84  13*38*56  FUTIL  6.18 


1! LENGTH  HRT  FP 
210 

310.5 

411.5 
512.99946 
6 '4.90839 
7 '6.99732 
8! 13.9946 
9*20.992 

10127.9893 

1 1 *34.9866 

12I4J.984 

13M8.9813 

14155.9786 

15162.9759 

16169.9733 

17176.9706 

18183.9679 

49190.9652 

20197.9625 

21 1 104.96 

22!  1 1 1 .957 

231 118.955 

241125.952 

251132.949 

261139.947 


OFFSET 

0 

1.79298 

3.41525 

5.08919 

6.77129 

8.1157 

11.36/4 

13.1284 

13.5 

13.5 

13.5 

13.5 

13.5 

13.5 

13.5 

13.5 

13.5 

13.5 

13.5 

13.334 

12.444 

10.7078 

8.07667 

4.51632 

0 


M IT. &BCD. POWER 

4-  9-84  10*13*59  FUTIL  6.18 


I ! SPEED  AND  PflMER  RESULTS  FOR  13.461  DESIGN  PROJECT* 
2! 

3! NON-PROPULS I ON  LOADS  <KN)t  135 
41PROPULSIVE  COEFFICIENT  *  0.77 
5! 

6! 


7'SPEED  CKT) 

Rl 

SHP 

TOTAL  KN 

9! 

1 

0.615661 

135.459 

10! 

2 

4.61783 

138.444 

11! 

3 

15.0496 

146.222 

12! 

4 

34.8408 

160.981 

13! 

5 

66.8572 

184.855 

14! 

6 

113.92 

219.95 

15! 

7 

178.818 

268.345 

16! 

8 

264.312 

332.097 

17! 

9 

373.141 

413.251 

18! 

to 

508.023 

513.833 

19! 

It 

671.663 

635.859 

20! 

12 

866.75 

781.335 

21! 

13 

1095.96 

952.256 

22! 

14 

1361.95 

1150.6! 

23! 

15 

1667.38 

1378.37 

24! 

16 

2014.9 

1637.51 

25! 

17 

2407.13 

1930 

26! 

18 

2846.71 

2257.79 

27! 

19 

3336.24 

2622.83 

28! 

20 

3878.34 

3027.08 

29! 

21 

4475.63 

3472.47 

30! 

22 

5130.68 

3960.95 

31! 

23 

5846. II 

4494. 44 

32! 

24 

6624.47 

5074.87 

33! 

25 

7468.37 

5704.17 

34! 

26 

8380.38 

6384.25 

35! 

27 

9363.07 

7117.04 

36! 

28 

10419 

7904.45 

37! 

29 

11550.7 

8748.37 

38! 

30 

12760.8 

9650.74 

39! 

31 

M051.8 

10613.4 

40! 

32 

15426.3 

11638.4 

41 ! 

33 

16886.8 

12727.5 

42! 

34 

18435.8 

13882.6 

43! 

35 

20075.9 

I5I05.6 

44! 

36 

21809.6 

16398.4 

45! 

37 

23639.5 

1 7762.9 

46! 

38 

25568 

19201 • 1 

47! 

39 

27597.7 

20714.6 

48! 

40 

29731.2 

22305.5 

-75- 


M!T. &BCD. BALnUT 
5-  6-84  10*54*3/  FUTIL  6.21 


1 ‘RESULTS  FROM  BALANCE  MODULE 
2! 

3!  GROUP 
4! 

5! 

6! 

/! 
til 
9! 

10! 

II! 


2 

3 

4 

5 

6 
7 

12!  Al 


WEIGHT 

LCG 

VCG 

860.354 

110 

15.5 

623.431 

120 

10 

90.814 

112 

15 

105 

60 

20.5 

243.00 

110 

0.5 

195.5 

/O 

22.5 

90 

90 

10 

22UJ.OO 

106.10 

13.002 

221.800 

76.889 

8.1414 

2430.9 

103.53 

13.378 

183. 64d 

00 

16 

327.043 

101.39 

15.5 

16!  MBT 
17! 

16!MARG1N  LEAD  (TONS)*  50  #VCGef>/2.  LCG*L«A/2 
INSTABILITY  LEAD  (TONS)*  171.809  CO*  6  FT 
20!ST ABILITY  LEAD  HEQUiRED  LCG*  67.1016 
21 ! 

22 ! NSC  (TONS)*  2623.55  MLCB*  103. 145 
23! 

24SSU9M  (TONS)*  2051.40  <LC3*  102.051 
Po'SUBMERGED  STABILITY  (RG)*  1.72272  FT 


MIT.  (WT.&BCD. POLY 
5-  7-84  10*24*02  FUT1L  0.21 


I! POLYGON  SUMMARY  OUTPUT  FILE 


2! 

3! 

4!  ITEM  LCG 

5! 

6!CREH  AND  EFFECTS  55 

/! SLUM'S  ON  COMP.  HATER  0 

INSANITARY  TANKS  85 

<3! LUBE  OIL  IN  SUMPS  00 

10'FIXED  CLEAN  FUEL  OIL  O 

1I!***«T<ITAL  FIXED  ITEMS 78.013 
I 2! PROVISIONS  AND  STORES  75 

!3!REVITALIZATI0N  02  88 

I4JT0HPED0S  IN  FHD  R<X)M  45 

155TORPEDOS  IN  AFT  Room  0 

16'TACTICAL  MISSILES  FHD  0 

1 7!TACTICAL  MISSILES  AFT  116 

I8!HHT  TANKS  30 

10! RESERVE  ELECTROLYTE  58 

20!TOTAL  FRESH  HATER  100 

21 ! RESERVE  LUBE  OIL  05 

22! FUEL  9ALLAST  Fhi)  26.5 

23!FUEL  BALLAST  APT  115 

24! VARIABLE  FUEL  OIL  45 

25! 

26! 

27! VARIABLE  BALLAST  TANK  OATA* 

28!  TANK  GR(XJP  LCG 

20!FoRHARD  TRIM  35 

30! AUXILIARY  /5 

31 ! AFTER  TRIM  105 


PTe**PP!I 


i 

t 


t 


£• 

/! 

N»N*I  '  1 

'&  } 

d! 

IF  (N.E0.3.OR.N.60.5)  GOTO  L(X)P  5 

a  r 

9! 

PRNT  DO  YOU  WANT  TO  SPECIFY  (S>  OH  CALCULATE  (C>  GROUP  IN]  WEIGHTS?  ( 

C  s 

10! 

read  ar  *  i 

|  / 

II! 

IF  (AH.F.O.  «S")  F(N)»I  I 

P  ' 

12! 

IF  (AR.EQ.-C-)  F(N)»0 

§ 

13! 

IF  (F(N).EO.O)  GOTO  Loop 

£ 

14! 

prnt  Enter  group  ini  weight  (tons).  „■ 

15! 

READ  NT 

16! 

H(N)=WT 

-18- 


MiT.ABCD.DKSIGN 

5-  7-04  13*50*00  FUTIL  o.2> 


I  !01M  H( 1 7) «P( 1 7 )  ,F(d).RS(2l  ).CIRC(2I 1  .LX(8>  .RX(d)  .XPdSOO)  .SA(2l )  <* 
2!>.LS<21).HP(4I ),KW(41).l)U(2n.X(29).Y(29> 

3 'HEAD  (ENTER  DESIGN  NAME*  )  ANAME 
4!«EICHT  CONTINUE 
5!N*l 
6! REPEAT 


I7!#L(X)P  CONTINUE 
ltt'UNTIUN.EQ.7) 

I9!REA0(maX.  OPERATING  DEPTH  IN  FEET*  »2*<AX 

20!P(  I  )s.2I55-M.606225E-4*2«AX 

21 !READ( RESERVE  BUOYANCY  FRACTION*  )P( (3) 

22 ! REaO( LEAO  FRACTION*  )P(Q) 

23!  RE  AD  <VAR  I  ABLE  LoaO  FRACTION*  )PUD 
24!READ<FREE  FL<X)D  FRACTION*  )P(15) 

25!PHNT  NUCLEAR  (II.  DIESEL  (2>.  OH  HANKEL  <3 >  PROPULSION? 
26! READ  PM 

27!PRNT  ATTACK  Cl)  OR  FBm  (2)? 

2d! READ  T 

2R!IF  (F(2>.E0.1)  GOTO  LABI 
30! COSUB  GHOUP2 

3 1 ! #LABI  IF  (F(6).EQ.n  GOTO  Al 
32! G0SU9  GR0UP6 
33! GOTO  A I 
34! •GH0UP2 

35!  READ  (SHAFT  HORSEPOWER*  )dHP 

3o!  IF  (PM.NE.U  GOTO  BATT 

37!  K*.03*SHP/I000*1.I3 

3d!  Nl»(K*l£5*SHP)/(LG(SHP)>**5 

3Q!  H(2>*Ni/2240 

40!RTNSUB 
4|!*BATT  CONTINUE 


VARTA  (2).  OR  TRIDENT  (3>? 


prnt  battery  type*  nicad  m 
read  bt 

RkADCToTAL  ENGINE  HORSEPOWER*  )tSP 
IF  ( PM.E0.2)  KE«.0096d 
IF  (PW.EQ.3)  KE*.(X)50d 
HE*KE*ESP 
WM».00645*SHP 

REaD(ToTal  kw-h  storage  for  sprint*  >ksp 

READCTOTAL  Kw-H  STORAGE  FOR  ENDUHANCE  HATE*  )KEMD 


57!  IF  (HT.EQ.3)  Cb*?5.4012 

5d!  BI=<ENQ/CE 

50!  IF  (Bl.CT.B)  B=Bl 

SI  PRNTEENTEH  DISCRETE  ADJUSTMENT  TO  GROUP  2  (0  OR  #  OF  TONS). 

62!  HEAD  DA 

63!  W(2)=WT*DA 

04!RTNSUB 

65 ! *GftOUP6 

66!  HEAD! NUMBER  OF  PERSONNEL*  )  NP 

6/!  rt(6)*2.3*NP 

6d!HTNSUB 

60!#AI 

70!  IF  (PM.EQ.1)  P(3)*.045 

71 j IF  t PM.E0.2.0R. PM.E0.3)  P<3)*.0I4 

72!?<4)=.06 

73! IF  (T.EO.I)  PC5)«.1I 
/4!1F  CT.EO. 2)  P(5)*.Ob5 
/5!IF  CT.EO.  1.AND.PM.E0. 1)  PC/)*. 03 
/6!IF  CT.EO. 1 .AND.PM.E0.2)  P(/>*.04 
/7! IF  (T.EO.I .AND.PM.E0.3)  PC  7)*. 05 
/d!lF  CT.EO. 2)  PC/)*. 14 

7Q»F«=(P<  l)*PCI)*PC0)>/Cl-PCm)4P(3)4PC5) 

BO! IF  (F(4) .EO.O. AND.FC 7) .EQ.O)  GOTO  CASEl 
tj  I ;  IF  CFC 4)  .£0.  I  .AND.  F(  7) .  EQ.O)  GOTO  CASE2 
B2!IF  CFC4).E0.0.AND.F<7).E0. I)  GOTO  CASE3 
B1!IF  CF(  4)  .£0.  I  .  AND.FC  7 )  .EQ.  I  )  GOT')  CASE4 
d4!#CASEl 

d5!  AI*C*CC2)^H(6)>/n-FH-P(4)-PC  /  )) 

d6!  W(4)*PC4)*A1 

d/!  WC/)=PC/)*AI 

dd!G0Tc>  TOTAL 

dO! #CASE2 

QO!  AI=CHC2)^N<4)^H(6))/< l-FH-PC 7)) 

O I !  rtC7)*PC7)*AI 
02! GOTO  TOTAL 
03I4CASE3 

04!  AI*CH(2)>N(6)«N( /))/CI-FH-PC4)) 

05!  NC4)=P(4)«A1 

06! GOTO  TOTAL 
07!#CASE4 

Od!  AI*(H(2)+H<4)*H(6)*W(  /))/(  l-F*) 

00! •TOTAL 
100!HC3)*PC3)*AI 
101 !HC5>*P<5)*AI 
102!HCd)*AI 
I03!HCO)»PCO)*AI 
IC)4!HC  IO)*AI+HCO) 

I05!H<  I2)*N(  IO)/C  I -PC  ID) 

I06!HC II >*P( II )*N( 12) 

I07!rt(l )*PC I )*NCI2> 

|Od!M( I3>*p( I3)*NC 12) 

IOO!4(  |4)=NC  12)*W(13) 

I I0!4CI5)*P<15)*H(I4) 

1 1  I !MC I6)*NC I5)4*< 14 ) 

1 1 2!PHNT  CALCULATES  COMPLETE ... 

I  l3!PRNf 
I I45PRNT 

1 15! <*  OUTPUT  SECTION 
1 1 0!  IF  CT.IG.  I)  <T*"4TlAr*V* 


IF  IT. HO. 2)  AT— BALLISTIC  MISSILE" 

IF  (PM.EQ.I)  &PM— NUCLEAR" 

IF  (PM.EQ.2)  APM— niESEL  ELECTH1C- 
IF  (PM.EO.3)  APM— HANKEL  ELECTRIC" 

1F  (DA.NE.O)  AM—  MOOIFIED  “ 

IF  (DA.EQ.O)  AM— 

IF  I3T.EQ.D  ABT— NICKEL  CADMIUM" 

IF  (BT.EQ.2)  ABT— LEAD  ACID  (VAHTA)- 
IF  I9T.EQ.3)  ABT— LEAD  ACID  (TRIDENT)" 

PHNT  HEIGHT  ESTIMATE  FOR  I&NAMEJ* 

PANT 

PRNT  TYPE*  (AT) 

PHNT  PROPULSION*  ( AM] (APM) 

IF  CPM.EO. 1 )  GOTO  SHAFT 
PHNT  BATTERY  TYPE*  I ABTJ 
PRNT  SPRINT  KH-H*  (KSP) 

PRNT  ENDURANCE  KH-H*  (KCNO) 

PRNT  BATTERY  HEIGHT*  IB) 

PRNT  ENGINE  HORSEPOWER*  (ESP) 

PHNT  ENGINE  HEIGHT*  (HE) 

PRNT  MOTOR  HEIGHT*  (HM) 

IF  (DA.NE.O)  PHNT  DISCRETE  ADDITION  TO  GROUP  2*  (DA) 

mshaft  continue 

PRNT  SHAFT  HORSEPOWER*  (SHP) 

PRNT  MAX.  OPERATING  DEPTH  (FT)*  tZMAXJ 
IF  (F(6).EQ.D  PRNT  NUMBER  OF  PERSONNEL*  INP) 

PRNT  LEAD  FRACTION*  f P (Q) 1 

PRNT  VARIABLE  LOAD  FRACTION*  [PHD) 

PRNT  RESERVE  BUOYANCY  FRACTION  *  (P(I3>) 

PRNT  FREE  FLOOD  FRACTION*  (P(l5>) 

PRNT  HIT  <RETURN>  TO  CONTINUE... 

READ  ACOMT 


150! PRNT  CATEGORY 

151 !PRNT 

I52M-1 

153! REPEAT 

HEIGHT 

154! PRNT  GROUP  (I) 
155! 1*1*1 

156!UNTIL  (I.EO.B) 

(HU)) 

1 5 /! PRNT  COND.  Ai 

(W(d) ] 

15B5PRNT  LEAD 

(H(Q)I 

IS^SPHHT  COND.  A 

(H ( 1 0) ] 

160! PRNT  VAR.  LOAD 

(NUI)) 

lollPRNT  NSC 

(H(12M 

162! PRNT  MBT 

(H( 13)1 

163SPRNT  SUB.  DISP. 

(HIM)) 

164! PRNT  FREE  FLO. 

(H( 15)) 

165JPRNT  ENV.  DISP. 
166! PRNT 

(H(lO)J 

PRNT  SATISFACTORY  (S)  OH  RECALCULATE  <R>? 

HEAD  A  DECK) 

IF  (ADECIO.EQ. "R")  GOTO  HEIGHT 

PHNT  PLEASF  HAIT  WHILE  OUTPUT  AND  PASS  FILES  AttF  HRITTEN.  THE  OUTPUT  For 
PHNT  THF  HbIGHT  CALCULATIONS  HILL  BE  IN  FILE  -MIT. HEIGHT-.  HHfcN  YOJ  HAVE 
PHNT  COMPLETED  THE  ENTIRE  DESIGN  PROGRAM,  "MIT. HEIGHT-  AND  THE  OUTPUT 
PRNT  FILES  FOR  S IBSEOUENT  MODULES  MAY  BE  OBTAINED  IN  HARD  COPY  BY  USING 
PHNT  THE  Command  'PRINTLP-FILFnamc— . 

OHtNH  2,-VlT.Ht  lGHT- 

*P  —  H 


T 


f!ASM?*SHP 

i!  4£S?*ESP 

'!AKSP«KSP 

HAKENO-KEND 

•AZMAX-ZMAX 

!ANP=Np 

U&BB-B 

«AWM«HM 

•ADA=DA 
!&LF=P(R) 
l'AVLF*P( II) 

!AR8F«P( 13) 
i!AFFF»PU5> 

!AW1=H(  I ) 

!!AH2=N(2> 

•AW3=W<3) 

!AW4*W<4> 

>!AW5*WC5> 

i!AW6=W<6) 

!4H7«W(7) 

!&H8«M(d) 
i!&WR*N(Q) 
l!  AW  IO*N{  10) 

!&H||»M(M) 

!!AW|2*W(  12) 

!AWI3»W( 1 3) 

,!4H|4»K(14) 
i!AWl5*W(  IS) 

>!AWI6*H(16) 

!AX«"HE1GHT  ESTIMATE  Foa  "♦4NA«E-*"*« 

! GOSUB  O 
tNRITEF  2.AP 
i!AX="TYP£*  "MT 
•CUSUB  O 

!! AX*** PROPULSION*  "♦AM*APM 
•IF  ( PM.EO. 1 )  GOTO  NUKE 
•COSU3  O 

,!AX*"BATTEHY  TYPE*  -♦ABT 
! GOSUB  O 

•AX--SPRINT  KW-H*  -♦AKSP 
ilCOSUB  0 

!AX»"ENDURAHCE  KW«H*  -♦AKEND 
l!  GOSUB  O 

•AX^-BaTTERY  HEIGHT*  -♦ABB 
IGOSUB  O 

;'4X*"ENGINE  HORSEPOWER*  -+AESP 
!G0SUB  O 

,!4X*"ENGINE  &  ASSOC.  EQUIP.  WEIGHT*  "♦AW£ 

•  !  GOSUB  O 

•AXe"MOToR  A  ASSOC.  EQUIP.  WEIGHT*  "♦AWM 
•GOSUB  O 

•IF  I0A.EQ.J)  GOTO  NUKE 

>!4X=*QiSCRETF  AQOITIon  TO  CROUP  2  WEIGHT*  -♦ADA 
•GOSUB  O 
•#NUKE  CONTINUE 

•AXS-4IAX.  OPERATING  DEPTH  IN  FEET*  -♦AZMAX 
•GOSUB  o 
JWH1TEF  2.AP 

»!AX*-IFAI1  FMAGTION*  "♦aLF 


J 


237! GOSUB  O 

P3BJAX— VARIABLE  LOAD  FRACTION*  -♦AVLF 
23Q5GOSUB  0 

240 ! 4 X— RESERVE  BUOYANCY  FRACTION*  "♦AHBF 
241 'GOSUB  O 

2421AX— FREE  FLOOD  FRACTION*  •♦AFFF 

243‘GOSUB  O 

244!HkITEF  2.AP 

2451NRITEF  2.AF 

?4o!  AX— CATEGORY  HEIGHT** 

24/1  GOSUB  O 
24a!wRIT£F  2.AP 
240'AX— GROUP  I  -♦4H1 

250!C0SJ3  O 

251  !AX— GROUP  2  -♦4H2 

2521GOSJB  O 

253!  AX— GROUP  3  *+l*3 

2541GOSU3  O 

255IAX— GROUP  4  "♦AH4 

256! GOSUB  O 

25714/— GROUP  5  -♦AHS 

25d! GOSUB  O 

25<J!  AX— GROUP  6  "♦AH6 

2601 GOSUB  O 
261  !6X— GROUP  7 
262JC0SUB  O 

263!AX— COND.  A I  "♦AHd 

264! GOSUB  O 

265  !AX— LEAD  -♦&H0 

2 66! GOSUB  0 

26/!  AX— COND.  A  -♦AH  10 

26d! GOSUB  0 

269! AX— VAH.  LOAD  -♦4H1I 

270!GOSUB  0 

271!  AX— NSC  -♦A**  1 2 

272! GOSUB  0 

2731AX— M2T  -♦AH  1 3 

2 /4 l GOSUB  0 

2751AX— SUB.  OISP.  "♦4H1H 

2/6!  GOSUB  O 

277!4X— FREE  FLO.  -♦AHlS 

2 /d!  G(1S‘JB  0 

279!AX— ENV.  OISP.  -+AHI6 

2 BO! GOSUB  0 
2B1!G0T0  DONE 
?tl2!#0  CONTINUE 
2d3!H«ITEF  2. AX 
2B41RTNSU9 
2tf5!»0ONE  CONTINUE 
2d6!OPENN  3.-NIT.PASS. " 

2«/!hKITEF  3.AH12 
2rtd!HRITEF  3. AH 15 
2691HR1TEF  3.AHlo 
2°0!HRirEF  3.ANANE 
2**  1 !  OPENH  4,-hIT.HALPaSS" 

292!HR'TEF  4.  AH1 
2Q31HHITSF  4,  AH? 

2N41HRITEF  4.  4H3 
295!hRITEF  4.  AH4 
JO'S* Hit ITFF  4.  *H«i 


[\tv 


M1T.&BCD.C3 

4-27-84  11*08i26  FUTIL  6.21 


I!  DU  RS<21>,ClRCC2l ),LXt8>,RX(«).XPOS<30>.SA(2l  J.LSC21  X# 

2!>*HP(41  ) .KM(4 1 )  aDU(  21  >  «X(  I32J.YC  132>.SX(21) 

3*#GE0M 

4'OPENH  1 .•MIT. PASS 1- 
SIREAOF  I «&M 
6!ONSO«H 
7!RFAUF  l.&N 
8!0FF«&W 
91HEADF  l.&H 
10!DENV*&* 

1 1 1READF  1 .ANAME 

1 25  P 1*3- 14159265 

I3140RANCH  CONTINUE 

14'READCEKTER  LENGTH  uF  ENTRANCE*  )LF 

I5!READ(ENTER  LENGTH  OF  HUN*  >LA 

16!REA0(ENTER  HULL  DIAMETER*  )D 

171REA0C ENTER  FORE9O0Y  PRISMATIC  COEFF.*  )CPF 

18!READ< ENTER  AFTERBODY  PRISMATIC  COEFF.*  )CPA 

19!L*LF»LA 

201LL-LF 

21!INTV»L/20 

22!R«0/2 

23!NF*71 .0477*CFF**5+65.8I U7*CPF**4-483. 3J37*CPF**3*587.4137*CPF**2<# 

24 • >-2d 1 . 7224*CPF*4Q .5876 

25!NA«379.6546*CPA**5-938.4708*C?A*^4+044.8853*CPA**3-471.0872*CPA**2<# 

26  5  >♦ 1 19. 1465 *CPA- 11.3454 

27!LS(U<J 

2tf!LS<2D»L 

29U-2 

30! REPEAT 

3t!LSIJ>*(J-l>MNTV 

32!XL»LS(J) 

33JGOSU3  RAO 
34!RSU)*RANS 
35!SA<J)*PI*RS«J)**2 
36!0U(J)*SA(J) 

37! 

38!UNTIL( J.E0.2I ) 

39! COSUB  INTO 
40 ! YEN  V»35*OEN  V 
41!DV*VENV-ANS 
42!LM8*DV/< P1*R**2) 

43!LL*LF'*LM3 
44!L0A*LF4LMB^LA 
45!lNTV*LOA/20 
46!LS(21  >«Loa 
47«LS< l>*0 
48! J*2 
49! REPEAT 

60!LSCJ)*(J-1 )*INTV 
51 ! XL*LS<  J> 

52!GOSUB  HAD 

53!RSCJ»*RANS 

54!CIHCCJ)*2*PI*HS(J) 

55!SA(J>*P1*HS«J>**2 
56!0U(  M=CIHC(J» 


-85- 


S7!J-J*1 

5d!UNTIL( J.E0.21 f 
59IG0SUB  INTO 
60!Sr«ANS 

6 1 ! VCYL*PI*LOA*R**2 
62JCP-VEMV/VCYL 
63!J*2 
64! REPEAT 

65!0U<JJ=LS(J)*SACJ) 

66! J=J*l 

678UNTJU  J.EQ-21 ) 

68JGOSU9  INTO 

60iLCffoUB«ANS/VEHV 

70!DR*.5*R 

71 !ON«ONSC*.8*PFF 

72!LO*DN-.005*0N 

73!Hl*0N».005*Df* 

M'PRMT  ON*  (ON) 

751PRNT  Uis  CLOJ 
76!  PR  NT  Hit  (Hi ) 

77!#LCBLP  CONTINUE 

7d!PRNT  ITERATING  ON  NSC  DRAFT  *  T=  HM)81 

/0!J*2 

dO! REPEAT 

81!IF<DR.C£.RSCJ)>  GOTO  SANE 
82!AC«(ACO$(DR/RS(J>))*(Pi/18l» 

d3!SXU)«SACJ)-(RS(J)*»2*AC-DH*(RSIJ>**2-DR**2)«*.5> 
84! GOTO  SKIP 
8S!#SANE  CONTINUE 
86!SXCJ)=SA(J) 

»7!#SKIP  CONTINUE 

88!DU(J)=SX(J) 

tt9!J-J*l 

90!UNTIL( J.EQ.21 ) 

QISCOSUB  INTO 
02!DN«ANS/35 

93!  IFCDN.GE.LO. AND. ON. (»£. HI )  GOTO  NON 
04JIFCDN.GT.HI)  GOTO  MINUS 
05!DR=DR^. I*DH 
96 'GOTO  LCBLP 
97!«MINUS  CONTINUE 
98!0R=0R-.I*DR 
99! GOTO  LCBLP 
1 00 ! #N()N  CONTINUE 
10 ! ! J“2 
102! REPEAT 

103!DU(J)*LS<J)*SX<J> 

104!  J*J*I 

I05!UNTIL( J.E0.2I ) 

I 06 ! GoSUB  INTO 
I07!LCBNSC*ANS/(0N«35) 

IOd!UX>UNTV-2)/5 
I09JLXC I )*.5 
110!LX(2)*1 
1 1 1 !LX(3)= 1 .5 
U2'LX(4)=2 

1 1 3! REPEAT  L«)P2--J*S,l.cj.GT.o) 

114'LXCJJaf J-4)«0X*2 
I !5!#LooP2  CONTINUE 
1 1 6! HfcPrAT  EXTRA:  1*1 .1 .(  I.GT.d) 


"8b- 


II  7!XL»LX<J) 

U81C0SUB  RAD 

M9!ftX(J)»RANS 

!20!#6XTflA 

»21!G0T0  GOJJT 

1 22! ARAD  CONTINUE 

1 23! IF( XL.LE.LF)  GOTO  FRONT 

124*  I FC  yUGT.LF.AND. XL.LE.  LL>  RANS*rt 

1 25! IF( XL.GT. LL)  GOTO  AFT 

I 26! GOTO  L(X)P4 

127!  CFRONT  CONTINUE 

1 28! XX*LF-XL 

1 29 ! RANS«H* ( ( 1 -< XX/LF) **NF ) •• ( »/NF ) ) 

130SC0TO  LOOP4 
131! #AFT  CONTINUE 
132!XX*Xl-U 

I33!RANS*R*C 1-(XX/LA)**NA) 

I34!#L(X)P4  RTNSUB 

I35!#INTC  <*S  UPSON'S  HULL  INTEGRATION 
I36!0D0*0 

1 37I2EPEAT  oD0l*J*2.2.<J.GT.20) 
l3d!oDO*onfH4*DU{  J) 

I30!#OOOL  CONTINUE 
140!£VEN=0 

141 ! REPEAT  EVENL*J*3*2.(J.GT.IR> 

1 42! EVENa£VEN^2«OU <  J ) 

143!4EV£NL  CONTINUE 

144!ANS*tINTV/31»lDU<l)*ory>*EVEN^DU(2l  )) 

l 45 'RTNSUB 

1 46 ! *GOUT 

147!L0*L0A/0 

148! PR NT 

149! PRNT 

1 50! PR NT  ENVELOPE  DISPLACEMENT  *  IDENV] 

15MPRNT  LENGTH  OVERALL*  ILOA) 

152! PRNT  OVERALL  PRISMATIC*  I  CP I 
153!PRNT  DIAMETER*  tDJ 
154!PRNT  LENGTH/DI AMEf ER  *  CLDI 
155! PRNT  MIDBO0Y  LENGTH*  {LN31 
lt>6!PRNT  FNO  PRISMATIC*  (CPF! 

157! PRNT  AFT  PRISMATIC*  (CPAl 

l50!PRNT  ENTRANCE*  ILFJ 

I59*PRNT  RUN*  CUI 

160* PRNT  LCB  SUBMERGED*  (LCBSUBl 

161!  PRNT  LCB  •  NSC  *  I LCBNSC1 

162! PRNT  DRAFT  •  NSC  *  IR^DRJ 

163!PRNT  SATISFACTORY  <S>  OR  RECALCULATE  Mil? 

164! READ  A0EC1D 

165!  IF  (AnECin.EO.«R">  GOTO  BRANCH 

166! PRNT  PLEASE  HAlT  WHILE  OUTPUT  AND  PASS  FILES  aRS  WRITTEN,  THE  OUTPUT 
167IPRNT  FILE  FOR  GEOMETRY  SUMMARY  DATA  WILL  5E  "MIT^EOMOUT'*.  A  FILE  OF 
J 68! PRNT  OFFSETS  WILL  BE  IN  -MIT.oCFoUT". 

I6R!0PENh  3.-MIT.OFFOUT" 

17O5A0—LENGTH  wRT  FP  OFFSET" 

1 7 l !GOSUB  W 
l72!LS(l)-0 
1  73!RS*11*0 
174!iX*LS(l) 

1  75!AC*RS(H 

1  /‘MAS***  " 


8 


—81— 


I 77!AQ=AX*AS+AC 
178!*RITEF  3.40 

1  79! REPEAT  L00P5*  J*l  .1  .U.GT.8) 
ltiO!AX-LX(J> 

t81!AC«RXCJ> 
lb2!AQ«4X*AS«-AC 
183! GOSUB  H 
184!J*J*I 

ld5!#Ux)P5  CONTINUE 

186! REPEAT  L<X)P6*J=2.I  .t J.CT.2I) 

ld7!4X=LS(J) 

188!4C*RS(  J) 

18<3!40»4X*4S*4C 
IQO'COSUB  * 

IQ!  !  ML00P6  CONTINUE 
1Q2!C0T0  FILE 
1Q3!*N  CONTINUE 
IQ4INRITEF  3.40 
1Q5! RTNSUB 
1Q6!4FILE  CONTINUE 
IQ/'OPENH  4. "MIT*  PASSHP" 

1Q8!4XcCP 

199‘NRITEF  4.4X 

200!4X*LOA 

20!!HRIT£F  4.4X 

202  !AXO 

203!NRITEF  4.4X 

204! AX*SF 

205!NRITEF  4.4X 

206!WRITEF  4.4NAME 

207JOPENH  5.-M1T.GEONOUT- 

208!  AX—GEOMETRY  OUTPUT  FOR  -♦4NAME*"*" 

20QINRITEF  5  .AX 

210!HRITEF  5.4S 

2JMWRITEF  5.4S 

2I2!4P»L0A 

2 1 3!  AX— LENGTH  OVERALL*  "♦AP 
2!4!HRITEF  5*AX 
21 5!4P*LMB 

2  t6!AX— MID-BUOY  LENGTH*  **4P 
217'HRITEF  5.4X 

218!4P«CP 

2  W!  AX— OVERALL  PRISMATIC*  "*AP 

220JNRITEF  5 .AX 

22l!AP*lO 

222!4X— LENGTH/DIAMETER*  "♦AP 
223!HRITEF  5, AX 
224! AP*LCBSUB 

225! AX— SUBMERGED  LCB  (AFT  ()F  FP)«  -♦AP 

226INHITEF  5, AX 

22/!rr*R*0R 

?28!AP*TT 

229!  AX— DRAFT  *  NORMAL  SURF.  C(>ND.*  "♦4P 

230!MRirEF  5, AX 

23MAP*LCBNSC 

2321AX— LCB  *  NORMAL  SURF.  COND.  *  -♦4P 

233!NRITEF  5*AX 

234INRITEF  5.4S 

235!AX— HSCA?  OF  INPUT  DATA*" 

236!WR1TEF  5.AX 


23/1MRUEF  5. AS 
23S1AP-DENV 

23<J!AX— ENVELOPE  DISPLACEMENT  —♦AP 

240!*RITEF  5. AX 

24I!AP«LF 

242! AX—  ENTRANCE*  -♦AP 
243!HRITEF  5. AX 
244 !  AP=LA 

245!AX— RUM  »  "♦AP 

246!HR1TEF  5. AX 
247!AP«D 

24tt! AX— DIAMETER*  •♦AP 
249!HRITEF  5.AX 
250! AP^CPF 

25 1! AX— F*D  PRISMATIC*  •♦AP 
2S2!HRITEF  5. AX 
253!AP*CPA 

254! AX— AFT  PRISMATIC*  "*SP 

255!NRITEF  5.AX 

256! OPEN*  6.-MIT.LC8" 

257!AX«LCBNSC 

25d!WRITEF  6. AX 

25<J!AX-LCBSUB 

2605MRITEF  6.AX 

261 !AX=R 

2625MR1TEF  6.AX 

263!AX»LOA 

264IMR1TEF  6. AX 

265!0PENN  7.-MIT.PASSPH1- 

266!  AX*D 

267!MRITEF  7. AX 

26tf!AX*LOA 

260!HRITEF  7, AX 

270!AX«LF 

27l!*RITEF  /.AX 

272! AX=LA 

273!NRITEF  7, AX 

274!AX*LCBSUB 

275!»»RITEF  7.AX 

2/6!DSUa«DENV-0FF 

277!AX«OSUB 

2/fl!WRiTEF  7, AX 

279!AX*rf>ENV 

2dC'»*RlTEF  /.AX 

2dl!AX*NF 

2&2!*RITEF  7.AX 

213!AX*NA 

2tt4!HRITEF  7. AX 

2d5!*0RAH 

2dd! X( I )*0 

20/!Y!U«0 

2ttd!>2 

2B9!REPcAT 

2«0!Xll»«LXCI-lI 

201  !YC  I)-RX(I-l) 

202!  !<■]♦! 

293!UNTIL< I.GT.9) 
2<*4*LFINT««LF-(1NTV+I  JJ/23 
2<J5!XL*INTV-1 
296!  >10 


REPEAT 
X<I)-XL 
GoSUB  RAD 
YUJ-RANS 
XL-XL+LFINT 
!•!♦! 

UMTIL(I.GT.3«> 

LINT*! LOA-LF)/ 101 

XL-LF+LINT 

1-32 

REPEAT 

X< I )*XL 

GliSUB  RAD 

YUI-RANS 

XL-XL+llNT 

1-1+1 

UNTIL! I.GT. 131 ) 

XI J32J-LOA 
YI1321-0 

OPEN*  B.-M1T. POINTS- 
J*1 

REPEAT 

AX-X(J) 

4Y-YCJ) 

NR1TEF  8.4X 
NRITEF  8.4Y 
J-J+l 

UNTIL! J.GT. 132) 

RUN  NEW  MIT.SPEED 
END 
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M  IT.  4BCD.  SPEED 
4-11-84  13*20*34  RJTIL  6.18 


IIDIM  HP(40),KH(40) 

2! READC ENTER  HOTEL  AND  COMBAT  SYSTEM  LOADS  IN  «H*  )HL 
3!READ(ENTER  PROPULSIVE  COEFFICIENT*  >PC 
4! OPENR  i »*M1T. PASSHP* 

5IREADF  I .4CP 
6! REAOF  l.&LOA 
7*READF  1.40 
8!READF  I  .ASF 
QSREADF  I, AN ARE 
10ICP-&CP 
I \ ! LOA«ALOA 
12!D»AD 
»3!SF-ASF 

1 4!  DA-!  .09065E-3*U>A*D*  1 1 .25 
15!J«1 
1 6! REPEAT 

1 7 ! RE - ( J  * 1 . 680%L0 A ) / 1 . 2  7  Q08E-5 
18!CF«.075/(LG(RE)-2J**2 

l0!CR«CF*CI.5*<D/L0A)**1.5^7*(n/L0A)**3)+.002*CCP-.6) 

20! CT*CF*CR*. 00025 

21  !HP(J)=(  1/PC)*.00872*J**3*<CT«SF*0A1 
22!  KH ( J ) *. 745  7*HP ( J  >*HL 
23!  J=J+I 

24! UNTIL* J.CT.40) 

25!PRNT 

26!PRNT 

27'PRNT 

28!PRNT  SPEED  AND  P0*»ER  OUTPUT  FOR  I&NAME ) 

20! PRNT 

30SPRNT  NON  PRoPULSlMi  LOADS  (K*0*  IHLJ 
3MPRNT  PROPULSIVE  COEFFICIENT  *  1PC1 
32!PRNT 

33! PRNT  HIT  <REHJRN»  TO  SEE  OUTPUT  VALUES... 

34JREA0  ACONT 
35! PR NT 

36'PRNT  SPEED  (KT)  SHP  TOTAL  K* 

37!PRNT 
38!J«1 
30 ‘REPEAT 

40! PRNT  CJ)  I HP( J>  J  IKM(J) ] 

41!J«>I 

42! UNTIL (J.GT.20) 

43! PRNT  HIT  <RETURN>  TO  CONTINUE... 

44 ‘READ  ACONT 

45! PRNT  SPEED  <KT)  SHP  TOTAL  KN 

46!PRNT 

4/!  REPEAT 

48! PRNT  IJ]  (HP(J)l  IKHU)1 

40! 

50!UNTIL(J.GT.40> 

5* !  PRNT  PLEASE  NAIT  HHILE  THE  OUTPUT  FILE  IS  CREATED.  OUTPUT  FOR  THIS 
52!PRHT  MODULE  HILL  BE  IN  FILE  •M1T.P0N£R«. 

53!o?ENH  2.-MIT.P0HEH- 

54!&X=*-SPEtO  AND  POHER  RESULTS  FOR  •♦ANAM£«-»- 

55!h«£TEF  2. AX 

56!A?ADl»* 


5/!4PA02«-  • 

58 SNRITEF  2.4PA0I 

50I4HL-HL 

60J4POPC 

4J14X"" NON- PROPULSION  LOADS  < KX »  t  ■♦AHL 
625NRITEF  2.4X 

63!4X»* PROPULSIVE  COEFFICIENT  i  -♦4PC 

64INRJTEF  2.4X 
65!NRITEF  2.4PA01 

66I4X«*SPEED  (KT)  SHP  TOTAL  KW« 

67'NRITEF  2.4PADI 

68JNRITEF  2,4X 

69JNRITEF  2.4PA0I 

70!  J«.t 

71 'REPEAT 

72!4J-J 

73!4X-HP(J) 

74!4HP«4X< 1.8) 

75!4Y-KH(J) 

73!4KH»4Y< I .8) 

77! 40UT**  •♦4J+APADI ♦"  "♦4HP+4PAD2VAKW 
78SHRITEF  2 ,4oUT 
79!  J*J+.l 

80!UNTIL( J.CT.40) 

81  .'RUN  NEW  MIT.  ENV 
82! END 
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MIT.&BCD.BALANCE 
4-27-84  10*47*27  FUTIL  6.21 


t  !DIM  M(  16 ) ,LC(  16)  .VC(  16) ,LM<  16)  ,VM(  16)  ,VY(16) 
2IOPENR  i.-MIT.BALPASS“ 

3!J*1 
4! REPEAT 
5IREADF  1 «&X 
6!N(J)«4X 


UNTIL(J.CT.l6) 

OPENR  2.-MIT.LC9- 
READF  2, AX 
LCBHSC-AX 
READF  2.&X 
LCBSUB*AX 
READF  2.&X 
R»AX 

READF  2, AX 
LOA*AX 

PRMT  THE  INPUT  DATA  FOR  THIS  MODULE  AND  THE  POLYGON  MODULE  MILL  BE 
PRNT  PROMPTED  IN  A  N£H  FORMAT.  WHEN  THE  PROGRAM  PROMPTS  FOR  AN  INPUT 
PRNT  VALUE.  THE  CURRENT  VALUE  HILL  BE  DISPLAYED  TO  THE  RIGHT  OF  THE 
PRNT  INPUT  PROMPT  MESSAGE.  THIS  WILL  APPEAR  AS  * 

PRNT 

PRNT  'INPUT  PROMPT*  «  CURRENT  VALUE' 

PRNT 

PRNT  THIS  WILL  FACILITATE  CHANGES  TO  INPUT  DATA  IF  MORE  THAN  ONE  ITER- 
PRNT  ATI ON  THROUGH  THE  MODULE(S)  IS  REQUIRED. 

45TART  CONTINUE 

LT-0 

VT»0 

J»l 

REPEAT 

PRNT 

PRNT  GROUP  [J]  LCG* 

REAOf VALUE?) LCIJ) 

PRNT  GROUP  CJ1  VCG* 

READ! VALUE?) VC(J) 

LM( J)»H( J)*LC(J) 

VM(J)«H(J)*VC(J) 

LT-LT+LMU) 

VT*VT>VM(J) 

J«J*I 

UNT1L(J.GT.  7) 

LC(d)*LT/N(8) 

VCC8)aVT/HC8> 

LM(8)*LT 

VM<8)*VT 

PRNT 

HEADCLCG  OF  VARIABLE  LOAD?)LC( 1 1 ' 

READ (VCG  OF  VARIABLE  LOAD?)VC<ll) 

PRNT 

LM( II >«M( 1 1 )*LC< II ) 

vm(  u>«N(in*vc(in 

PRNT  TOTAL  lead  IS  IN(*IJ  Tons. 

PRNT 

HE  AO  (AMOUNT  OF  LEAD  FOR  MARGIN  IN  T«>NS7>ML 
SL*N(<7)-ML 


i. 


f 
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57 •  LSL*  C HU  2 )  *LCBNSC-LM(  8  > -LM ( 1 1  )-ML* <  LOA/2 )  >/SL 
58 !  LC  (9  )«<SL*LSL*ML*  (LOA/2  ))/*<<>> 

591PRNT  REQUIRED  STABILITY  LEAD  LOG  IS  [LSLJ. 

60! HEAD! IS  THIS  FEASIBLE  <Y/N>?  )4R 
61 1  lF(4R.E0.*N->  GOTO  START 
62!LM(9)»W<9)*LC<9) 

63!VN(9>*SL*6*ML*R 

64!VC(9)*VM(9)/W(9> 

65!LM(10>*LM(8>*LM(9) 

66!VMU0)*VM(8>«-VMt9> 

67!  LC(  10)»LM(  10)/H(  10) 

68!VCU0)*VM<  10)/H(  10) 

69!LN()2)*l«( 10>*LMC 1 1 ) 

70! VH( 12)“VM{ 10)+VM< 1 1 ) 

71 ! VC( 1 3 )  *R 
721 VM ( I 3) *H*H (13) 

73!LC( 13)*(W( |4)*LCBSU8-LM( 12))/H( 13) 

74'PRNT 

75IPRNT  REQUIRED  MBT  LCG  IS  ILC( 13) ). 

76! READ! IS  THIS  FEASIBLE  <Y/N>?  )4R 
77! IF<4R.E0."N«)  GOTO  START 
7d!U<l3>«H<I3>*LCn3> 

79!LM<|4)*LM( 12)*L¥(I3) 

80!LC( I4)*LM( 14)/HC 14) 

8l!VM(»4)«VM(l2)*VM(l3> 

82! VC( 14>=VM(14)/H(14) 

83!BG*R-VCCU) 

84 ! PR NT 

85IPRNT  SUBMERGED  STABILITY  <8G>  IS  (BGI  FT. 

86! READ( IS  THIS  ACCEPTABLE  <Y/N>?  >4R 
87!  IF<4R.EC.  -N-)  GOTO  START 

88! PUNT  THE  OUTPUT  FOR  THIS  MODULE  HILL  BE  IN  FILE  -MIT. BALOUT-. 

89!J«I 

90! REPEAT 

91!VY(J)»VC(J)-R 

92!J*JU 

93!UNTIL(J.GT.14) 

94!oPENH  3,-MlT.BALOUT- 
95!4X«-RESULTS  FROM  BALANCE  MODULE- 
96!HRITEF  3,4X 
9/!4Pl  —  - 

98!4P2*" 

99!4P3*-  - 

1001HRITEF  3.4P1 

101 !4X*-GR0UP  HEIGHT  LCG  VCC- 
1021HRITEF  3.4X 
I03IHRITEF  3.4P1 

104! J«| 

I05JREPEAT 
106!4J*J 
I07!4X*H( J) 

I 08!4H*4X (1*0) 

10<3!4X«LCUJ 
110!4LC*4X<1 ,6) 

1 1  l!4X«VC(J> 

1  I21&VOAXU  .6) 

M3J40UT*-  -*4J+- 
IMIHR1TEF  3«40UT 
1 131 J»J*l 
1 16! UNTiL< J.GT. O 


'♦AH«-AP2«-4LC-*4P3*4VC 
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1  l7!J*d 

118!AG="  A1  « 

I 19JG0SUB  10 
I20!J*0 

121'AG-"  LEAD  - 
122SG0SUB  10 
I23!J*10 
124! AG*"  A  * 

125'COSUB  In 
126! J*l I 
1 27 !AG«"  V  LD  " 

1 281GOSU8  in 
I2Q!J*I3 
1 3Q!AG*M  KBT  - 
131‘GOSUB  10 
I 32! ANSC*W( 12) 

I  33!ASUB*H<14) 

134«&LCBN*LCBNSC 
1 35! ALCBS»LCBSUB 
136JABG-BG 
1371MRITEF  3.AP1 
138JASL-SL 
13<?!4kL*ML 

1 40! AoUT**MARGlN  LEAO  (TONS)*  "♦AML+"  *VCG*D/2.  LCG*Loa/2" 
1 4 1 !HR ITEF  3.40UT 

142!4nirr="STA3ILITY  LEAD  (TONS)*  "♦ASL-*"  6VCG«  6  FT" 
143!HRITEF  3 ,4 OUT 
M45  4LSL-LSL 

I45!40UT«"STABILITY  LEAD  REQUIRED  LCG*  "♦ALSL 
I46IWRITEF  3.A0UT 
I4/JHRITEF  3.4PI 

1 48!AOUT*"NSC  (TONS)  i  "♦4NSC*"  #LCB*  "♦ALCBN 
14QJHRITEF  3.40UT 
150IWRITEF  3.4P1 

151 !40UT*"SU3M  (TONS)*  "♦4SUB*"  «LCB*  "♦ALCBS 
1 52!  HR  ITEF  3.40UT 
153'HRITEF  3.API 

l 54 ! AOUT*" SUBMERGED  STABILITY  (BG) *  "♦ABC**  FT" 

155! HR ITEF  3.A0UT 
15610PENH  4  »"MIT»  PASSPOLY" 

157!AX=LC8SUB 
158!HRITEF  4, AX 
159!AX*LC( 11 ) 

160!HRITEF  4 ,4X 
161 !AX=LC( 13) 

I62!HRITEF  4. AX 
l 63! GnT n  PLOT 
164!# in  continue 

I  65! AX*H( J ) 

1 66!AH*AX( I ,8) 

1  6  /!  AX»LC  ( J) 

1 68 !4LC*AX ( I  ,6) 

16R!AX«VC(J> 

1 70! AVC*AX( 1  ,6 ) 

171  !A<HfT»AG*4H*A?2'*ALC+AP3*AVC 
1  /2'HRITEF  3,A(HJT 
173IRTNSUB 
1 74!#PL0T 

1 751SELECT  MODE  MODEL 

i/6! select  layer  3 


I  771 J»l 
I  781  REPEAT 
J7R1AJ-J 

1 80! INSERT  POINT* XILC(J) )Y(VY<J> J<CR> 

181 !XP«LC(J)*1 

1 82! INSERT  TEXT  t&JJ  HGT  1 .5* XtXP)YtVY(J) J<CR> 

183!J*J*1 

1 84! UNTIL {J.GT*  7) 

185! INSERT  POINT* XfLC(9>  JYCVYO)  ]<CR> 

186! XP*LC(9)+I 

187! INSERT  TEXT  'LEAD*  HGT  l.5*X[XPIYtVY(9) )<CR> 
1 80 ! INSERT  P0INT*X1LC( 13) 1 Y[VY< 13) J<CR> 

1 86! XP«LC( 1 3 )♦ I 

100! INSERT  TEXT  'MBT'  HOT  1 .5* XI XP1YI VY< 13> )<CR> 
I Q I ! RUN  NEW  HIT. POLY 


MIT.AHCO.POLY 

5-  7-34  16*14 *02  FUT1L  6.21 


IDiM  W<  16>,XVU0).YVU0>.XPU0).YP<l0).XTU(».YTt  10)  ,M  (o)  ,Y(d) 

I  OPENS  1 ."MIT.BALPASS- 

!J»1 

!  REPEAT 
IHEADF  |«ftX 
!W(J)=AX 

I  UNTIL<  J.GT. 16) 

!VLOaW< II) 

I‘4BT*W(I3) 

!DSUB=W( |4) 

! OPENS  2.MMlT. PASSPOLY" 

IREADF  2. AX 
!LSUB*&X 
IREADF  2. ax 
!LVLO*AX 
IREADF  2, AX 
ILMBT«&X 

•START  CONTINUE 
8AL=VL0 

PRNT  TOTAL  VARIABLE  LllAD  IS  IVLOl  TONS. 
iPrtNT 

READ< HEIGHT  OF  CREh  aNO  EFFECTS7)CE 

REAi)(LCG7)LCE 

:bal=bal-ce 

cosub  stat 

READ( HEIGHT  of  sLBm'S  OR  COMPENSATING  hater?) mc 

REa»(LCG?)LMC 

BaL«BAL-MC 

cosub  stat 

RE A!)C WEIGHT  OF  SANITARY  TANKS  ANO  WaTER?)ST 

REA0(LCG7)LST 

9aL*SAL-ST 

COSUB  STAT 

REAOtwFIGKT  OF  LUBE  OIL  IN  SUMPS?>lO 

NFAI)tl.CG?)LLO 

BAL*BAL-LO 

COSUB  STAT 

REAOCWEIGHT  OF  FIXED  CLFAN  F.  o.  <IE.  SHIELD  TANK>?)CF() 

RFAD(LCC?)LCFo 

BAL*RAL-CFo 

GOS'JB  STAT 

READC WEIGHT  OF  PROVISIONS  ANO  STORES? IPS 

RbAO(LCG?>LPS 

bal*hal-ps 

GOSUB  STAT 

READ! WEIGHT  OF  REVITALIZATION  OXYGEN 7) 02 
READ < LOG? >L02 
9aL»BAL-02 
GOSUB  STAT 

rfao(hf:ight  of  torhetos  jn  forward  hoom?>tf 
•lEAOCI  CG7ILTF 
BAL*3AL-Tr 
GiteUB  bTAT 

REAIXrtEIGHT  oF  TU'PEODS  I »i  AFT  RO0M7ITA 
RfcAiHI  CG7ILTA 


Cu.Vu'a. 


V-.  *.»•  w 


!9*L=BaL-Ta 
scosua  ST AT 

!READ(W£1G4T  OF  TACI1CAL  MISSILES  FoRMARD?)MF 

!ftEAD(LCG?)LMF 

!9AL*dAL-MF 

!G()SU6  STAT 

JkEADtwEICHT  OF  TACTICAL  MISSILES  AFT?>MA 

!REaD(LCG?)LMA 

!RAL»3AL-MA 

!COSU9  STAT 

!TM»TF-*TA-*Mr*«A 

:LTM=<TF*LTF*rA«LTA*MF*LVF^MA*LMM/TM 
1  HE ADt HEIGHT  OF  WATER  FoH  WRT  TANKS? )WRT 
!  RE  AD(LCG?>LWR7 
!BAL*8AL-wftT 
1GOSU9  STAT 

!3EAU(rtEICHT  OF  RESERVE  ELECTROLYTE?) EL 

!REAU<LCG?)LEL 

!3AL*8AL-EL 

1GOSUB  STAT 

!H£A0<  WEIGHT  OF  FRESH  WATER  <PoTT  ABLE. FEED .BATTERY  RES£HV£>?)Fw 

!HEAD(LCG?)LFw 

3AL«BAL-Fw 

1COSJB  STAT 

HEADtwEIGHT  OF  RESERVE  LUBE  4  HYDRAULIC  01L?>«Lo 
!REAO(LCG?)LRLO 
!  BAL*BAL-RLo 
i  COSUB  STAT 

IPHNT  THE  LAST  VARIABLE  LoaO  ITEM  TO  BE  CONSIDERED  IS  FUEL  OIL.  TO  ENSURE 
PRNT  A  MARGIN  of  SAFETY  AT  THE  BASE  OF  THE  POLYGON*  A  VARIABLE  FUEL  OIL 
PRNT  CVFO)  TANK  IS  AUTOMATICALLY  SPECIFIED.  THE  REQUIRED  VFO  CAPACITY  IS 
PRNT  COMPUTED  USING  YOUR  INPUT  OF  THE  FRACTION  0=  TOTAL  FUEL  TO 
PRNT  PLA^E  INTO  THE  VFo  TANK.  THE  MINIMUM  SUGGESTED  FRACTION  IS 
PRNT  .23  (THE  RATIO  DIFFERENTIAL  oF  FUEL  AND  WATER  DENSITIES). 

PRNT  IT  IS  HIGHLY  RECOMMENDED  THAT  AN  ADDITIONAL  MARGIN  ClF  AT 
PRNT  LEAST  .05  BE  SPECIFIED  TO  Alio*  A  SAFETY  MARGIN  AT  THE  BASE 
[PRNT  OF  THE  Polygon. 

PRNT 

Rfc AO (FRACTION  oF  FUEL  To  de  PLACED  IN  VFo?)  VFR 

7Fo*9AL 

VFo=VFH*TF() 

FBT»7Fo-YF() 

FL«CE*MC«-3T/3*LO*CF(> 

LfL=(CE*LC£-*MC*LMC*ST*LST*Lo»LLO»CFo*LCFo)/FL 
PRNT  FBT  CAPACITY  (TONS)*  IFBTI 
PRNT  VFO  CAPACITY  (TONS)*  IVFoJ 
PRNT 

HEAD( AMOUNT  OF  FBT  <TONS>  IN  FORWARD  7ANKS7) F9TF 

FBTA*FBT-FBTF 

PRNf 

REAf)(FwD  FBT  LCG?)LFB7F 

RE  A.)  (AFT  FBT  LCG7>LF3TA 

LFBT=( FHTF*LFHTF^FBTA*LFBTA)/FBT 

READ (VFO  LCG?)LVFo 

PRNT 

PRNT 

PRNT 

PRNf 

PRNT  VARIABLE  LOAD  RECAP* 

PRNT  IffcM  WEIGHT  LCD 


mrr* , 

. ■/.  , 
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I  I  /!  PRNT 

I  Id! PRNT  CHEW  A  EFFECTS  ICE) 

II9SPHNT  SUM'S  CMC) 

120JPRNT  SANITARY  TANKS  1ST) 

12HPRNT  SUMP  LUBE  OIL  tLOJ 

1 22! PRNT  CLEAN  F.O.  (FIXED)  ICFO) 

123iPRNT  ***TUTAL  FIXED***  IFLI 

1 24! PRNT  PROV.  A  STORES  IPS) 

125'PRNT  HEVITALI2ATION  o?  (02) 

I  26!  PRNT  TOHPEDOS  FwD  R(X)M  [Tr) 

127'PRMT  TOHPEDOS  aFT  ROOM  ITA) 

1 2d! PRNT  TACT-  MISSILES  PND  IMF) 

1 2^! PRNT  TACT.  MISSILES  AFT  IMA) 

130! PRNT  wht  TANKS  IHRT) 

131 'PRNT  RESERVE  ELECTROLYTE  I  EL) 

1 32! PRNT  TOTAL  FRESH  HATER  (FW) 

1 33! PRNT  RESERVE  LUBE  OIL  [RLoJ 

134 'PRNT  FUEL  BALLAST  F*0  IFBTF) 

1 35! PRNT  FUEL  BALLAST  '*=T  I F3TA) 

1 36 1PRNT  VARIABLE  FUF  .  IVrO) 


I3/!REA0(ARE  THESE  VaLUcS  ACCEPTABU  <Y/N>? 
I38!IF(AR.EO.*NO">  GOTO  START 
1 70! PRNT 

l 40 • AFL»LSUB-LFL 

141!APS«LSUB-LPS 

142«A02=LSUB-Lo2 

143!ATM*LSUB-LTM 

I  44 • AWRT3LSUB-LHRT 

I4S!AEL*LSUB-LEL 

146!AFH*LSUB-LFH 

|47!ARLO*LSUB-LHLO 

14d«AFBT*LSUB-LFBT 

14Q!AVFO*LSUB-LVFO 

150!aTF*LSUB-LTF 

151 ! ATA*LSUB-LTA 

152!AMF*LSUB-LMF 

I53!AMA»LSUB-LMA 

I54!aFBTF-LSUB-LFBTF 

I  55!AFBTA»LSUB-LFBTA 

I 56 ! MFL*FL*AFL 

157!MPS*PS*APS 

158!M02*02*A02 

15Q!MTM«TM»ATM 

I 60!MHRT*HRT*AHRT 

161 !MEL*EL*AEL 

162!MFHaFH*AFH 

I63!MRU»*RL0*ARL0 

1 64 ! MFBT*FBT*AFBT 

165!  MVFo*VF()*AVFO 

1 66! UTF*TF*ATF 

I67!MTA«TA*ATA 

1 6dl MMF*MF*AMF 

l6*i!MMA*MA*AMA 

1 70! «FBTFaFBTF*AFBTF 

I  /l!«FtJTA*F9TA*AF3TA 

172!CVL*.QV375 

I /3 ! CVH* I • 0063 /5 

i 74! Goto  nsTm 

I /5!#STaT  CONTINUE 

I  76! PHNT 


ILCE) 

ILMC) 

ILSTJ 

I  LLo) 

ILCFoJ 

ILFL) 

ILPS) 

(L02) 

ILTFI 

ILTA) 

ILMF) 

I  LMA) 

(LhrT) 

(LED 

ILFW) 

ILRLO) 

I LFBTF ) 
ILFBTAJ 
ILVFO) 
)4R 


1  7/JPHNT  BALANCE  C»F  VARIABLE  LOAD  ACCOUNT*  tBALl 

17d!PRNT 

1  70'HTNSUB 

I0O!#NSTM  CONTINUE 

tdI!LDSUB«W(14)-*<IO) 

182!#HEAV1  CONTINUE 
1835CHM3/35 

1 8  4  •  H 1  *FI>TM«- .  S*PS+£L+H  RT+ .  S*  RLO* .  5  *  FW+CH  *FBT 

1 85!  *H I -HFL+MTM*.  5**PS*VEL*MWRT*.5»NRLo«-.  !>*MFW*C**MFBT 

Id6! VBHl -LDSUB-MBT-H 1 

1 87! AM8T»LSUB-LN3T 

I  titsi  m*1BT»MBT*AMBT 

I  dO !  4VBH1  *-MHBT-«Hl 

|O0! XV(I )*NVBH1*CVL 

101 !YV( 1 )*VBH1*CVL- 

10?!#HEAV2  CONTINUE 

1 03  •  H2*FUTM*Ps*£L*WRT*RUHFW^C«*FBT 

104!MH2*MFL-Hfr«^NPS^MEUMNRT-H*RL(HMrW>O(*MFBT 

105! V8H2=L0SUB-MBT-H2 

106! NVBH2*-MMBT-NH2 

|07!XV(2J*HVBH2*CVL 

I Od ! YV ( 2 ) * V8H2*CVL 

100!*LIT6I 

200 ! L 1  *FI>. 75*PS«*HT^.  75*RL<H.  75*Fh+FBT*VFO 

201 ! NL 1  *MFL+.  75*MPS*M»«HT*.  75*#RLO*.  75*NFh«*FBT*MVFO 

202!V8LI»LDSUB-NBT-L: 

203! MVBL 1 *- MMBT-N  L 1 
204'XV(3)»MVBL!*CVH 
205!YV(3»*VBL1*CVH 
206!#L1TE2 

20 7 !  L2*FL*  NRT+RLO*  VFlHCN  *FBT 

208 !  ML2«NFUM*«RT*MaL<HM  VFO*CW  *  «F3T 

200 !  VB  L2»L0bU  BH4BT-L2 

2 1 0 !  NV  BU  »-NMBT-NL2 

211! XV (4) =NVBL2*CVH 

212!YV(4)-VBL2»CVH 

213!#HFI 

2M!HF1«FL<»MF*TF*.  75«PS*EL*Hin>.  75*RLO*.75*FW«-VFo*FBTA-K:h»FBTF 

2I5!MHF1»MFL«*TF^MMF*.  /5*MPS«-MEL*MNRT«-.  /5*MRU»«  /5*nFH*mVFo*MFSTA*CW*VFBTF 

2 1 6!  VBHF1  »U)SUB-MBT-HF1 

21 7! NVBHF I *-MMBT-NHF1 

2I8!XV(5)«NVBHFI 

210'YVC5)*VBHFI 

220!  #...-2 

?21!Hr2*FL*MF4TF*.5*P3*EL*Hirr*.25*RL<H.5*FN^VFo»FBTA*CN«FBTF 

222!MHF2*NFL*MTF*W1F^»75#NPS-*NEL*nhRT'*.25*MRLo^.5*NFW*MVFo*MFBTa+CW*mFBTF 

223 ! VBHF2*LDSU  B-MBT-HF2 

224  !*vBHF2«-N»«BT-MHF2 

22b!XV(6)«4VBHF2 

226 ! YV (6 ) *VBHF2 

227!*HA 

223!HA*FL-fMA^TA^*b*PS^HHT^.75«HUH-.5*FH+VFlHFBTF^CH*FBTA 

22Q!MHA«»FL*‘«A+‘*TA+.5*HPS^NWRT^.75«NRU)^.5*ilFN*NVRHIiFBTF-K:4«MFBTA 

230 ! VBHA*L0SU8-MBT-HA 

231 !NVBHA“-NMBI-NHA 

232 ! XV( 7 ) =NVBHA 

233! YV ( / ) *VBHA 

234!#N 

23b!N=H2 

236!«M*NH2 


23  /! V9N»L0SUB-MBT-N 

238UVBN— MMBT-MN 

239!XV<8>*MVBN 

240!YV(d)-V9N 

241 !#CONOM  CONTINUE 

242!CM*H2*CH*FBT+FBT 

2431 MCM*MH2-C**MF8T*KFBT 

2  44 ! V8CM-LDSUB-MBT-CM 

245UVBCM—MMBT-MCM 

246!XV<9J*MVBCM 

247!YV<9)«V8CM 

243! <#•****•*•*****  SET  UP  LOCATION  OF  TRIM  A  AUX  TANKS  *******•••••*** 
249!PHNT 

250! PRNT  YOUR  ACTUAL  DESIGN  MAY  HAVE  AS  MANY  TRIM  AND  AUX  TANKS  AS  YOU  DEEM 
251 !PRNT  APPROPRIATE  FOR  YOUR  PARTICULAR  ARRANGEMENT  SCHEME.  TO  SIMPLIFY  THE 
252! PRNT  CALCULATION  PROCESS  FOR  THE  EQUILIBRIUM  POLYGON,  YOU  APE  ASKED  To 
253! PRNT  ENrER  AGGREGATE  CAPACITY  AND  LCG  FOR  THE  FORWARD  TRIM  GROUP.  THE 
254! PRNT  AFTER  TRIM  GROUP.  AND  THE  AUXILIARY  GROUP  TANKS. 

255!PRNT 
256! GOTO  SIZE 
257!#L<X)P 
25U!EXECV 

25R!READ(LCG  OF  FORWARD  TRIM  TANK  CROUP7)LFT 

260!READCLCG  OF  AFTER  THU  TANK  GROUP ?)LAT 

26UREA0CLCG  OF  AUX  TANK  GHOUP7JLAUX 

262JAFT-LSUB-LFT 

263!AAT*LSUB-LAT 

264JAAUX-LSU9-LAUX 

265!RTNSUB 

266!#S1ZE  CONTINUE 

267! J«1 

263! REPEAT 

269! IFC XV( J).GT.5000)GUT<>  LARGE 
2 /0! IF(XV (J).LT. -5000) GOTO  LARGE 
27I!IF(YV(J).GT.500)GDTO  LARGE 
272!J*J*» 

273!UNTIL( J.GT.9) 

274 'GOTO  SMALL 
2/5ULARGE 

276!DEL  PAR  MIT.HULLOUT 
277! EXIT  PAHT  FILE  MIT.HULLOUT 
2  7d ! XS* 1 000 
2  79!YS*100 

280! ACT  PAH  MIT.POLYOUT 

281! ACT  ORA  POLYGON  FORM  MIT.LaHGEPOLY  DRAW  PoLY 
282! SEL  CPL  LEFT 

2d3!0EF  VIE  GRAPH*  XII Y6.X0Y0.X22YI 7  <CR> 

284! GOTO  GRAPH 

285MSMALL  CONTINUE 

286! DEL  PAR  MIT.HULLOUT 

287! EX IT  PART  FILE  MIT.HULLOUT 

280!XS*500 

289! YS*50 

290! DEL  PAR  MIT.POLYOUT 
29 I! ACT  PAR  MIT.POLYOUT 

292! ACT  ORA  POLYGON  FORM  MIT.SMALLPOLY  DRAW  POLYP 
293ISEL  CPL  LEFT 

294 !DfcF  VIE  GRAPH*  XII Y6.X0Y0.XP2YI /  <CR> 

295!#GHAPH  CONTINUE 
296 !J* I 
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297! REPEAT 

29B!XP(J)*XV(J)/XS 

2Q9!YP(J)*YV(J)/YS 

300! INSERT  POINT*  XI  XP( J) IYI YP( J)  1  <CR> 

30| !  XT<  J )  -XP(  J  )♦.  2 
302 !  YT  X  J )  -YP ( J )  -.  1 
303!J«J*I 
304!UNTIL(J.GT.9) 

305!  INS  TEX  Hi  HC  .2*  XIXTU  I IYIYTC1  )  I  <CH> 

306! INS  TEX  H2  HC  .2*  X(XT(2) JYl YT(2) 1  <CH> 

307! INS  TEX  Ll  HG  .2*  XI  XT(3)  J  YI  YT  (3>  1  <CR> 

30BJ1NS  TEX  L2  HC  .2*  XI XT(4) J YI YT(4 ) J  <CR> 

309! INS  TEX  HFl  HG  .2*  XIXTC5) JYl YT(5) 1  <Cfl> 

310! INS  TEX  HF2  HG  .2*  XI XT(6 ) JYC YT( 6 ) 1  <CR> 

311 ! INS  TEX  Ha  HC  .2*  XI XT( 7)1 YI YT( 7) ]  <CR> 

312! INS  TEX  N  HG  .2*  XI XT<d>-.5 J YC YT(B> I  <CR> 

313! INS  TEX  M  HG  .2*  XlXT(9) IYIYTC9) J  <CR> 

3 MJGOSUB  LOOP 
3I5JGOTO  PLTNK 
316!#CAP  CONTINUE 
3I7!EXECV 

3 IBiHE AO (CAPACITY  <TONS>  OF  FORWARD  TRIM  CROUP?) CFT 

3 19! READ (CAPACITY  <TONS>  OF  AUXILIARY  GROUP?)CAUX 

320!REA0( CAPACITY  <TONb>  ()F  AFTER  TRIM  GR0JP7JCAT 

32l!RTNSUB 

322! #P LINK  CONTINUE 

323JG0SUB  CAP 

324!M(l)*0 

325!M(2)-AFT*CFT 

326 !M (3 )*M( 2 )+AAUX*CAUX 

327!M(4>«M(3)*AAT*CAT 

328!M(5)«M(4)-aFT*CFT 

329 ! M( 6) -M(5) — AAUX*CAUX 

330!M(7)*0 

33 1 S YC 1 )»0 

332!Y(2)-CFT 

333!YU)=CAUX*CFT 

334! Y(4)«Y(3)^CAT 

335!Y(5)*Y(4|-CFT 

336!<C6)*Y(5)-CAUX 

337!Y(/>*0 

33(5  !J*1 

339JREPEAT 

340!M( J)«M( J)/XS 

341 !Y( J)*Y( JJ/YS 

342!J«J«-l 

343!UNTIL(J.GT.7> 

344!  INS  LIN  TAG-1  *  XIM(|)  JYIY(in.XlM(2)iY[Y(2>J,  <CR> 

345!INS  LIN  TAG-2*  XIM(2) I YC Y(2)1.XI MC3) IYI Y(3> ) •  <CR> 

346!  INS  LIN  TAG-3*  X(N<3MY!  Y(3>  1  .XfM(4)JYI  Y«4)  J,  <CR> 

347! INS  LIN  TAG-4*  XIMt4)]Y(Y<4)) «X(M(5) IYI Y(5) 1 »  <CR> 

348! INS  LIN  TAG-5*  XIM(5) IYI Y(5> ) «XIM(6) IYI Y(6) 1 ,  <CR> 

349! INS  LIN  TAG-4*  XtM(ft) ] YI Y(6) J .XIM( 7) IYI Y( 7 ) J .  <CR> 

35u!EXECV 

351 !READ( IS  THIS  POLYGON  SATISFACTORY  <Y/N>  ?  )AR 

3$2!lF(4H.EO.*Y")  GOTO  FINAL 

353! DEL  ENT*  TAG  1 .TAG  2.TAG  3.TAC  4, TAG  5. TAG  6  <CR> 

354!£XECV 

355!  READ  (CHANGE  CAPACITIES  ONLY  <l>.  OR  LCG'S  6.  CAPACITIES  <2>  7  >CHC 
356! IF(CHC.E0.2)  GOTO  BOTH 


357! 1F<CHC.£Q. I )  GOTO  PLTNK 

35S!#8oTH  CONTINUE 

359JGOSUB  L(X)P 

3oO!GnT()  PLTNK 

36 1 !#FINAL  CONTINUE 

362JOPENW  3  »*M IT. PASS.  P OLYDAT 

363J&X-CE 

3645GOSUB  0 

365!4X-LCE 

366! COSUB  O 

367!&A«NC 

368'GOSUB  O 

369!&X=LMC 

3/OJGOSUB  0 

371 !4X*ST 

3/2JGOSUB  0 

373!4X*LST 

374‘GOSUB  O 

37514X-LO 

376IGOSUB  0 

3  7/!4X*ao 

37d!GoSUB  O 

3  /9!&X*CF0 

3dO!G(>SUB  O 

381 !4X»LCFO 

382IGOSUB  0 

3d3!4X«FL 

3d4!G0SUB  0 

385!4X*LFL 

386IGOSU3  0 

3«7!&X*PS 

3881GOSUB  O 

3d9!AX»LPS 

3QOJGOSUB  O 

39|'AX*02 

392JGOSUB  O 

393!AX*l02 

394! COSUB  0 

395!&X«TF 

396‘GClSUB  0 

39/!4X«LTF 

398!GOSUB  0 

399J4X-TA 

400!G0S’JB  0 

401 !4X*LTA 

402JGOSUB  0 

403J*X»I*F 

404JGCSU1  O 

405!4X*U»F 

406!GOSUB  O 

407!AX*‘«A 

4Gd!G<)3Uo  O 

409:4<*l«a 

4*0!  Gt -SUB  0 

4 1 | ! AX*HRT 

4J2JGOSUB  () 

413!AX*LHHT 
41 4  !G(lSUB  O 
4I5!AX»EL 
416 !G0SU9  O 


4l7!4X*L£L 
41ti!GOSUB  « 

4|Q!4X*FH 

42U1GOSUB  O 

421UX-LFW 
422! GoSUB  0 
423!4X*HLO 
424!G0SUB  O 
425!4X*LHLO 
426JCOSUB  a 
427!&X-FBTF 
428*.GOSUB  0 
429!&X«LFaiF 
430! GOSUB  O 
43» !AX«FBTA 
432!G<)SUB  O 
433!4X*LFBTA 
434! COSUB  » 

435!4X*VFO 
436! GOSUB  n 
437!&X*LVRI 
4 3d! COSUB  I) 

43<J!4X»CFT 
440!GoSUB  » 

441 !&X»LFT 
4421GOSUB  () 

443!iX«CA0X 
444 ! GOSUB  C) 

445!&X*LAUX 
4461GOSU3  » 

447! AX-CAT 

44d! COSUB  O 

44Q!&X«LAT 

450JGOSUB  O 

451! PLOT  DOT  SCALE  -5 

452 ‘EXIT  PART  F 

453! GOTO  END 

454 !#0 

455INR1TEF  3*AX 

450!RTN3UB 

457!«END 

458! RUN  NEW  MlT.PoLYFILfc 
459! END 
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MIT. A  BCD. POLYFILE 
5-  7-84  JQ«22*03  FUTIL  6.21 


IJOPENR  I  ."MIT.  PASS.  POLYDAT" 

2!OP£NH  2 ."M  IT.  OUT . POLY- 

31  PRNT  OUTPUT  FOR  THE  POLYGON  MODULE  KILL  BE  IN 
4 ! AT**PoLYGoN  SUMMARY  OUTPUT  FILE" 

51KRITEF  2. AT 
6!AP«" 

7!KRITEF  2.AP 
81HHITEF  2.AP 

Q!AT="  ITEM  LCG 

IO'kRITEF  2.  AT 
M'HRITEF  2.AP 
1 2 !  RE  ADF  I  .AX 
I3!READF  I. AY 

I4!AT*"CREH  AND  EFFECTS  -*AY*" 

I5IKRITEF  2. AT 
|6!G0SUB  Rf) 

I 7!AT»"SLBM'S  OR  COMP.  WATER  "♦AY*" 

I8!HRITEF  2. AT 
IQ!GOSUB  RD 

20!AT*"SAN1TAHY  TANKS  *•♦**♦" 

2HKR1TEF  2. AT 
225GOSUB  RO 

23!AT*-LUBE  OIL  IN  SUMPS  "♦AY-** 

24‘KRITEF  2. AT 
25! COSUB  RD 

26!AT—FIXED  CLEAN  FUEL  (UL  -♦AY+- 

27'KRITEF  2 .AT 
2U!G05UB  Rl) 

29*AT*"****ToTAL  FIXED  ITEMS****  "♦AY** 
30'KRITEF  2. AT 
31 ! COSUB  RD 

32! AT="PROV ISIoNS  AND  STORES  “♦AY*" 

33JKRITEF  2. AT 
34'GOSUB  RO 

35!AT***HEVITALIZAT10N  02  "♦AY*" 

36!KRITEF  2. AT 
37JCOSUB  RD 

38!AT*"T0RPED0S  IN  FHO  R(X)M  -*AY*" 

39!NfiITEF  2, AT 
40IGOSUB  Rl) 

41!AT*"TORPEOOS  IN  AFT  ROOM  “♦AY*- 

42!KRITEF  2, AT 
43! COSUB  RD 

44! AT* “TACTICAL  MISSILES  FHO  -♦AY*" 

45!HRITEF  2 .AT 
40!G0SUB  RD 

4/!AT*"TACTICAL  MISSILES  AFT  -♦AY** 

48'KHITEF  2. AT 
49!GOSU8  RO 

50! aT*"HRT  TANKS  "♦AY*" 

•jIMHlTFr  2, AT 
52!GOSUB  HD 

a3!AT*-RESERV6  ELECTROLYTE  "*AY*" 

54!HR1TEF  2.AT 
55IGOSUB  RO 

56! AT=-ToTAL  FRESH  WATKP 


'MIT. OUT. POLY'. 

HEIGHT" 

"♦AX 

"♦AX 

"♦AX 

"♦AX 

"♦AX 

"♦AX 

"♦AX 

"♦AX 

"♦AX 

"♦AX 

"♦AX 

"♦AX 

"♦AX 

"♦AX 


"♦AY*' 


"♦AX 


5 /!HRITEF  2. AT 
56SG0SUB  HD 

59! AT*" RESERVE  LUBE  OIL  "♦AY*"  "♦A* 

60! WH1TEF  2  .AT 
6HG0SUBRD 

62! AT* "FUEL  BALLAST  PHD  “♦AY*"  “♦** 

63!*RIT£F  2. AT 
64SGO5U0  RD 

65!AT*-FUEL  BALLAST  AFT  “♦AY**  '•♦AX 

66!*&ITEF  2. AT 
67'GOSUB  RD 

od! AT—  VARIABLE  FUEL  OIL  “♦AY*"  "+AX 

69'WRITEF  2 .Af 
70!G0SUB  HI) 

/HHHITEF  2.AP 
721KHITEF  2.AP 

/3! AT*"VARlABLE  BALLAST  TANK  DATA*" 

/4!HHITEF  2. AT 

75!AT“"  TANK  GROUP  LOG  CAPACITY  (TONS* 

76!*RITEF  2, AT 

77 !  at  ""Forward  trim  “♦ay*-  '•♦ax 

7d!HRlTEF  2, AT 
79! COSUB  HO 

dO! AT*" AUXILIARY  “♦AY*"  "*AX 

61 JNRITEF  2.AT 
62 JGOSUB  RO 

83!AT»"AFTER  TRIM  “♦AY*"  '•♦AX 

64!WRITEF  2. AT 

65! GOTO  END 

86!#RD  CONTINUE 

67! REAOF  I. AX 

6d!REA0F  1.AY 

691RTNSU3 

90!#END 

9i !  end 


